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INTRODUCTION 

Flexible  conduits  are  usually  designed  to  guard  against  five  primary 
modes  of  failure.   As  illustrated  in  Figure  1,  these  are:   (1)  wall 
crushing,  which  occurs  when  compressive  stresses  due  only  to  circumferential 
thrusts  exceed  the  yield  stress,  (2)  separation  of  seams,  when  the  thrust 
exceeds  the  seam  strength,  (3)  initiation  of  buckling  from  an  essentially 
elastic  state  of  stress,  (4)  inelastic  buckling,  and  (5)  excessive 
deflection,  or  flattening,  due  to  plastic  yielding  under  combined 
compressive  and  bending  stresses. 

It  is  comparatively  easy  to  guard  against  the  first  two  modes  of 
failure,  because  the  magnitude  of  the  circumferential  thrust  is  relatively 
less  sensitive  to  soil  properties  or  to  pipe  stiffness.   However,  bending 
stresses,  which  play  a  special  role  in  failure  modes  (3),  (4)  and  (5),  are 
very  sensitive  to  both  the  above  parameters,  and  simplified  procedures  to 
account  for  their  effects  -  for  example,  by  controlling  the  crown 
deflection  -  are  not  generally  applicable.   Good  progress  has  been  made  in 
recent  years  in  treating  culvert  problems  on  a  rational  basis  (Leonards 
and  Roy,  1976;  Wenzel  and  Parmellee,  1976;  Katona  et  al,  1976;  Duncan,  1977), 
Significant  advances  have  been  made,  and  more  can  be  expected  in  the 
future. 

This  report  deals  with  the  performance  limits  of  buried  flexible 
conduits,  with  special  emphasis  on  the  buckling  failure  modes.   A 
literature  review  was  made  to  trace  the  development  of  the  present 
state-of-the-art  for  predicting  buckling  of  flexible  conduits.   The 
assumptions  and  models  used  in  the  theoretical  analyses  are  discussed 
and  evaluated,  and  their  applicability  to  the  problem  at  hand  is 
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assessed.   Comparisons  of  predictions  of  buckling  pressures  by  the 
various  theories  for  some  typical  culvert  geometries  and  soil  strengths 
are  made  to  show  the  relative  effects  of  tlie  assumptions  usrd    In  i  lie 
respective  prediction  models.   Relevant  experimental  test  results  arc- 
reviewed  and  assessed  in  terras  of  the  test  conditions,  insofar  as  these 
reflect  the  manner  in  which  culverts  behave  in  situ.   Finally,  the 
problems  associated  with  predicting  failure  conditions  for  buried 
conduits  are  reviewed,  highlighting  those  areas  in  need  of  further 
research. 

DEFINITION  OF  BUCKLING  TERMS 
In  accordance  with  currently  accepted  definitions:  buckling  is  the 
development  of  a  kink,  wrinkle,  bulge,  crease,  or  other  sharp  change  in 
original  shape;  the  critical  load  (thrust,  stress,  or  external  pressure) 
is  the  load  at  which  two  alternative  load-deflection  relationships  are 
mathematically  possible,  i.e.  the  theoretical  load  at  which  buckling  is 
initiated  under  ideal  conditions;  the  buckling  load  is  the  load  at  which 
a  compressed  element,  or  member,  buckles  in  service  or  in  a  loading  test; 
an  instability  is  reached  when  the  capacity  to  resist  additional  load  is 
exhausted  and  continuing  deformation  results  in  a  decrease  in  load- 
carrying  capacity.   Because  real  structural  members  have  minor  geometric 
imperfections,  and  loading  is  never  perfectly  symmetrical,  buckling 
loads  are  generally  somewhat  lower  than  the  critical  loads.   Experiments 
have  shown  that  this  difference  is  especially  marked  for  unsupported 
thin  shells  subjected  to  hydrostatic  compression;  accordingly,  it  is 
important  to  keep  in  mind  that  buckling  theories  inherently  tend  to 
overpredict  the  buckling  load. 


4. 


Buckling  of  flexible  conduits  can  manifest  itself  in  a  number  of 
ways.   Elastic  buckling  is  initiated  at  stress  levels  below  yield. 
Buckling  in  a  high  mode  involves  a  large  number  of  small  "waves" 
distributed  around  the  circumference.   Buckling  in  a  low  mode  also 
involves  the  entire  circumference  but  is  characterized  by  a  small  number 
of  larger  waves.   Local  buckling  is  a  crimp  or  crease  involving  a  small 
percentage  of  the  conduit  circumference.   Because  the  mechanism  of 
failure  is  different,  the  several  types  of  elastic  buckling  occur  at 
different  critical  pressures;  however,  due  to  geometric  imperfections  and 
the  presence  of  residual  stresses,  even  in  the  case  of  buckling  in 
theoretically  high  modes  crimping  of  the  conduit  wall  may  be  visible  only 
at  one  or  two  locations.   For  this  reason  it  is  difficult  to  distinguish 
between  the  modes  of  elastic  buckling  solely  by  visual  inspection. 

Inelastic  buckling  is  initiated  after  portions  of  the  conduit  wall 
have  undergone  plastic  yielding  and  generally  occurs  in  a  low  buckling 
mode.   A  snap-through  buckle  involves  a  sudden  reversal  of  curvature 

or  inversion   in  the  conduit  wall  and  could  result  in  an  instability  if 
a  sufficient  fraction  of  the  circumference  is  involved.   On  the  other 
hand,  although  undesirable  from  a  design  standpoint,  local  buckling  may 
not  immediately  involve  instability. 

1.   BUCKLING  THEORIES 
The  analyses  presented  in  this  section  treat  buckling  of  thin-walled 
cylindrical  shells  or  rings.   The  buckled  shape  is  taken  to  be  constant 
along  the  axis  of  the  cylinder,  and  thus  no  longitudinal  buckling  modes 
resulting  from  axial  forces  are  considered.   Although  some  solutions  for 


longitudinal  buckling  are  available  (Almroth,  1962;  Reynolds,  1962; 
Weingarten,  1962),  to  be  consistent  with  the  fact  that  analyses  for 
predicting  performance  of  culverts  generally  assume  plane  strain 
conditions,  the  buckling  theories  reviewed  herein  are  concerned  with 
the  two-dimensional  problem  in  which  sections  transverse  to  the  axis 
of  the  pipe  are  assumed  to  behave  identically. 

A.   Elastic  Buckling  in  High  Modes 
Timoshenko  and  Gere  (1961) 

An  early  model  used  to  approximate  buckling  of  flexible  conduits 
buried  in  soil  was  a  circular  ring  subjected  to  a  uniform  hydrostatic 
external  pressure.   The  solution  for  the  critical  external  normal 
pressure,  p  ,  was  derived  by  Timoshenko  and  Gere  (1961)  and  given  as. 


R 


where,    E  =   Young's   modulus 

1  =   moment   of   inertia   of    ring  wall 

R  =    radius 

n  =  buckling  mode  =  i\R/ i,    where  £  =  half  wave  length  of  the 
buckled  shape 
When  solved  in  terms  of  the  critical  compressive  stress,  f  ,  within  the 
ring,  eqn  (1)  becomes,  in  the  lower  bound  (n  =  2) 


f  =  ^-^  (2) 

^        R^t  ^^ 


t  =    ring   thickness 
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The  above  solution  resulted  from  consideration  of  the  equilibrium  of  a 
deformed  ring  element  in  which  only  circumferential  stresses  were 
considered.   Nonlinear  deflection  terms  were  also  neglected.   Application 
of  this  solution  to  large  variations  in  ring  dimensions  showed  that 
incongruous  results  were  obtained  for  some  selected  dimensions  because 
of  the  linearization  simplification. 

Brockenbrough  (1964) 

Brockenb rough  (1964)  proposed  the  use  of  eqn  (2)  for  determining  the 
critical  stress  of  deeply  buried,  flexible  culverts.   The  depths  for 
which  the  analysis  was  deemed  applicable  were  to  exceed  those  at  which 
live  loads  on  the  soil  surface  had  a  significant  effect  on  the  culvert. 
From  field  experience  a  depth  of  ten  feet  was  chosen  as  a  limiting  value, 
independent  of  the  culvert  span.   Brockenbrough  believed  that  for  the 
above  conditions  the  compressive  stress  in  a  culvert  could  be  calculated 
from  the  expression  for  simple  hoop  compression,  as  suggested  by  Wtiite 
and  Layer  (1960), 


P  R 
f  =  -^  (3) 


f  =  circumferential  compressive  stress 

p  =  vertical  overburden  pressure  at  crown  =  Y  z 

where,  y  -   total  unit  weight  of  soil 

z  =  height  of  soil  cover  above  the  crown 

The  compressive  stress  calculated  by  eqn  (3)  could  then  be  compared  to 
the  critical  stress  obtained  by  eqn  (2)  to  determine  the  maximum  allowable 
height  of  soil  cover  for  a  given  conduit  stiffness  and  soil  density. 
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Approximating  the  uniform  external  radial  pressure  by  the  nominal  over- 
burden pressure  (yz)  was  thought  to  be  conservative  when  considering 
flexible  conduits,  because  soil  arching  was  apt  to  reduce  the  normal 
pressure  reaching  the  conduit  to  a  value  lower  than  yz . 

In  developing  his  suggested  maximum  fill  heights  above  buried 
conduits  of  given  thickness  and  corrugation  profile,  Brockenbrough 
reasoned  that  the  weight  of  the  overlying  soil  would  deflect  or  flatten 
the  culvert  and  cause   redistribution  of  soil  pressure  to  a  uniform 
value,  but  he  did  not  indicate  how  much  deflection  was  necessary  to 
develop  the  full  redistribution.   He  also  realized  that  Timoshenko  and 
Gere's  solution  did  not  consider  flattening  at  the  crown  and  invert, 
which  would  decrease  the  buckling  resistance.   Brockenbrough  thought  that 
the  use  of  a  "conservative"  equivalent  uniform  pressure  (p  =  yz)  would 
more  than  offset  the  neglect  of  initial  deflections.   In  addition,  he 
believed  that  localized  yielding  of  the  outer  fibers  of  the  culvert  wall 
due  to  combined  bending  and  compressive  stresses  would  probably  not 
cause  failure.   Thus,  when  the  ring  stress  calculated  in  eqn  (3)  reached 
the  critical  stress  formulated  by  Timoshenko  and  Gere,  the  conduit 
would  fail  by  elastic  buckling. 

Valentine  (1964) 

Valentine    (1964)   proposed  an  analysis   to   estimate   the 
critical   buckling   load  on   shallow    (less    than   one   diameter  of   soil   cover 
at   the  crown)    flexible  conduits   subjected  to   live   loads, 
assuming  that   a   portion   of   the   top   of   the   conduit   acted  as   a   two-hinged 
curved   bar.      With  this  assumption,    the   formula    for  arch   buckling   under  a 
uniform  normal   pressure    (Roark,    1954)    could  be  used   to   determine   the 


critical  pressure  acting  on  an  appropriate  section  of  the  conduit,  which 
could  be  represented  by  an  equivalent  arch, 


P^  " 


EI 


^    R^ 


2 

a 


(4) 


p  =  uniform  critical  pressure  acting  on  the  arch  section 
2a  =  central  angle  forming  the  arch,  in  radians,  0  <  a  <  y 


The  critical  pressure  calculated  by  eqn  (4)  was  correlated  to  an  equivalent 
surface  load  for  design  considerations.   The  distribution  of  load  on  the 
conduit,  reflected  by  a,  was  found  to  be  a  function  of  the  depth  of 
burial  and  conduit  radius  for  a  given  surface  loading.   Some  suggested 
values  of  a  were  back-calculated  from  a  series  of  field  tests  on  shallow 
buried  conduits  up  to  54  inches  in  diameter,  surrounded  by  dense  sand  and 
subjected  to  the  AASHO  H-20  design  loading;  when  the  depth  of  soil  cover 
was  one  radius,  a  equalled  0.62  and  0.69,  respectively,  for  54  in.  and 
24  in.  diameter  conduits.   To  be  conservative,  Valentine  suggested  the 
use  of  a  safety  factor  equal  to  two  with  the  analysis.   He  concluded  that 
buckling  in  conduits  with  shallow  burial  was  critical  for  typical  conduits 
if  the  diameter  was  greater  than  42  inches.   However,  no  limiting  soil 
depth  was  suggested  to  differentiate  shallow  from  deep  burial,  no  general 
expression  for  calculating  a  was  given,  and  the  validity  of  the 
equivalent  uniform  pressure  distribution  used  by  Valentine  is  questionable. 
Nevertheless,  this  marked  one  of  the  first  attempts  to  differentiate 
between  the  behavior  of  shallow  and  deep  conduits. 

Bodner  (1958) 

Bodner  (1958)  extended  the  von  Mises   theory  of  shell  buckling,  which 
calculated  the  critical  buckling  pressure  of  a  long  cylinder 
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subjected   to   uniform  normal   pressure,    to    consider  uniform  constant-direction 
boundary  pressures.      As   in   the   von   Mises   solution,    biaxial    stresses  were 
considered   and   inextensional   buckling,    in  which   the  neutral   axis   of   the 
shell  wall   remains   constant    in   length,    was   assumed.      However,    Bodner's 
solution  allowed   the  direction   of   the   boundary   pressures   to    remain 
constant,    while   the  normal   pressures    in   the  von   Mises'    case  were    forced   to 
change  direction   in   order   to    remain   perpendicular   to   the   deflecting   shell. 
The   resulting   expressions    for  the   critical   buckling   pressures   included   a 
factor,   J,   not    included   in   the   ring  buckling   expressions. 


1-v 


p     =  — ^ r—      {von  Mises}  (5a) 

""        R-^(l-V^) 


p     =  -^ TT      {Bodner}  (5b) 

"     r\i-v2) 


V  =   Poisson's   ratio 

From  equations    (5)    Bodner  concluded  that,    in   general,    a   uniform  constant- 
direction   pressure   is   less   critical   than   a   uniform  normal   pressure. 

Armenakas   and   Herrmann    (1963) 

Development    of   buckling   analyses    continued  with  the  work  of   Armenakas 
and  Herrmann    (1963),    who   derived  a  more   general   solution   that   considered 
buckling  of  cylinders  of  finite  length  under  uniform  constant- 
direction  and  unifoLTii  normal   pressures.      The   effect   of   the   buckling   mode 
on   the   critical  pressure  was   also   accounted   for.      In   the   analysis,    the 
final   deformed  configuration  was   attained  after  the   shell   passed   through 
an   intermediate  equilibrium   state  of    initial    stress    from  an  original 
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unstressed   condition.      In   the  equilibrium  equations,    initial   displacements 

accompanying   the   initial    stress   state   were   included,    and   it   was    found   that 

neglect   of    these   displacements    in   previous   analyses   produced   some   error. 

However,    the   initial    displacements   could   be  neglected  without    error    in    the 

equations   of  motion   determining   the   buckled  configuration,    and   this   allowed 

a   closed   form  solution    for  critical   normal   and  constant-direction   buckling 

pressures.      Armenakas  and   Herrmann   found   that    for   finite   length  cylinders 

with  high  buckling  modes,    the   critical   constant-direction   buckling  pressure 

was   essentially   equal   to   the   critical   normal   buckling  pressure.       In 

cylinders  with  large   length   to    radius   ratios    (—  >   50)   and   relatively   low 

R 

buckling  modes,    the   critical   constant-direction   buckling  pressure  was 
greater   than   the   critical  normal   buckling  pressure,    as    given   by, 


Pc  =   ^*  -^h-3  <^) 

(1-V   )R 

2 

K''^  =  for  constant-direction   pressure 

»/ 1' 

K*  =  (n"  -  1)  for  normal  pressure 
n  =  buckling  mode 

Equation  (6)  reduces  to  eqn  (5a)  when  the  K'^  expression  for  normal  pressure 

is  used  and  n  =  2.   For  thin  shells  (—  <<  1),  eqn  (6)  also  reduces  to 

R 

eqn  (5b)  when  n  =  2. 

The  general  buckling  equation  for  —  <  50  showed  that  the  critical 

K 

buckling  pressure   for  a   finite  length  cylinder  was   greater   than    that    for 
a   cylinder  of   infinite   length.      Thus,    the   use  of   eqn    (6)   will   give 
conservative   estimates    for   the   critical   buckling  pressures,    provided 
the   loading  and   deformations   are   strictly   two-dimensional. 
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Donnell  (1956) 

Donnell  investigated  the  lack  of  agreement  between  experimentally 
measured  external  buckling  pressures  of  cylindrical  shells  and  theoret- 
ically calculated  critical  pressures.   Experimental  buckling  pressures 
were  often  much  less  than  the  theoretically  predicted  values,  and  the 
observed  buckling  modes  were  much  lower  than  the  predicted  modes. 
Although  these  discrepancies  had  been  attributed  to  geometric  and 
material  imperfections,  no  attempt  had  previously  been  made  to  relate 
the  decrease  in  buckling  resistance  to  the  degree  of  imperfection. 
Donnell  considered  imperfection  effects  on  shell  buckling  by  modifying 
the  extended  von  Mises  theory  considering  finite  shell  lengths  to  in- 
clude an  unevenness  factor  in  the  load-deformation  equations.   The 
unevenness  factor  represented  the  expected  geometric  imperfection  in  an 
unloaded  shell,  which  is  dependent  on  the  shell  material  and  Llic  inaiui- 
facturing  process.   The  imperfection  was  represented  by  an  initial  strain, 
and  the  critical  pressure  was  then  calculated  as  a  function  of  the  shell 
geometry,  yield  strength,  elastic  modulus,  and  unevenness  factor,  U. 

The  results  showed  that  the  normalized  critical  buckling  pressure 


varied  with  geometry  and  imperfection  level  as  follows: 


(1)  the  buckling  pressure  decreases  as  the  shell  length  increases, 
for  any  given  imevenness  factor,  U. 

(2)  As  the  flexibility  (— )  of  a  shell  increases,  the  influence  of 
the  imperfection  level  on  the  critical  buckling  pressure  becomes 
less  pronoxinced. 
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Figure  2    shows  a   set   of    typical    results    from   Donuell's   solutions    for 

long   tubes,   >   9.      The    important    influence  of    the   yield   stress   on    the 

buckling  pressure   is   clearly    indicated.      Tests   performed  at    the   David 
Taylor  Model   Basin   showed   back-calculated  uneveness    factors    to    range 
typically    from  U  =    0.0001    to    0.0005    in   smooth  steel    tubes   with  —  between 
75  and  250. 


Cheney    (1963) 

Cheney  (1963)  analyzed  the  stability  of  a  circular  ring  under  plane 
stress  conditions,  subject  to  circumferential  support  by  elastic  springs 
and  acted  upon  by  a  uniform  external  pressure  on  the  ring  wall  of 
constant  magnitude  and  direction.   Initial  strain  prior  to  buckling  was 
neglected.   Cheney  attempted  to  derive  a  more  general  solution  by 
retaining  nonlinear  energy  terms  throughout  the  analysis.   His  results 
showed  that  circumferential  support  of  the  ring  forced  it  to  buckle  in 
higher  modes  and  at  higher  pressures  than  the  unsupported  case.   Although 
Cheney's  representation  of  support  by  elastic  springs  was  an  over- 
simplification, it  opened  new  avenues  of  analysis  that  could  model  more 
closely  the  effect  of  soil  restraint  on  the  critical  buckJ.ing  prt  ssure. 
In  the  solution, a  uniform  external  pressure  distribution  was  assumed  to 
develop  around  the  cylinder  after  very  small  or  negligible  ring  deformation. 
This  was  necessary  to  obtain  a  closed  form  solution,  as  given  by  eqn  (7). 


P^  =  -7  (n   -  1)  +  -^ (7) 

R  n  -1 


_3 
k  =  coefficient  of  elastic  reaction,  or  'subgrade  modulus'  (FL   ) 
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FIGURE  2    CRITICAL  PRESSURE  AS  A  FUNCTION  OF 
GEOMETRY  AND  IMPERFECTION  LEVEL 
(after  DONNELL,  1956) 
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Forrestal  and  Herrmann  (1965) 

Forrestal  and  Herrmann  (1965)  extended  Cheney's  analysis  to  consider 
cylindrical  shells  subjected  to  a  uniform,  constant-direction,  external 
pressure  on  the  shell  wall  and  supported  by  a  continuous,  elastic  medium 
rather  than  discrete  springs.   The  theory  of  elastic  buckling  in  the 
presence  of  initial  stress  and  displacement  was  used  in  the  general  solution 
of  the  complicated  boundary  value  problem.   At  the  critical  pressure,  p  , 
a  deformed  equilibrium  position  was  obtained,  which  induced  surface 
tractions  on  the  shell  and  caused  changes  in  the  initial  uniform  pressure 
distribution.   These  surface  tractions  were  related  to  the  shell  displace- 
ments, and  equations  developed  by  Armenakas  and  Herrmann  (1963)  were  used 
to  determine  the  expression  for  the  critical  pressure  causing  buckling 
which,  for  the  fully  bonded  case,  is  given  by. 


2   2  y^^'-V 

n  (n  -l)(3-4v  )/3  +  ^ — 

■^      3EI/r''(1-v^) 


2n(l-v  )  +  (l-2v  ) 
s         s 


2 
n  t 


(8) 


(n  -l)(3-4v  )  +  2n(l-2v  )  4-  1 
s  s 


p      =   critical,    uniform  hydrostatic   pressure   for   buckling  of   imsupported 

cylinders,    as   expressed   by   eqn    (5a) 
E     =   Young's  modulus   of   supporting  medium 
V     =   Poisson's   ratio   of   supporting  medium 
A  solution  was  also   obtained   for  no    interface    friction,    as   given   by. 


p  2    _  E  /(1+v   ) 

c  _   n   -1  s  s 

Pq  ^  3El/R-^(l-v^) 


(1-2V   )(n+l)+n 


(9) 
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The  minimum  (or  lower  bound)  buckling  pressure  for  an  elasticall   s  pported 
cylinder  with  given  support  and  stiffness  properties  is  obtained  by  ^^ub- 

stitution  of  buckling  mode  numbers  into  the  above  equations  until  the 

Pc 
minimum  —  ratio  is  calculated,  in  effect,  minimizing  eqns  (8)  and  (9) 

o 

with   respect   to  n. 

Equation    (8)    shows   that   the  critical   buckling   load   increases   only 
15%  as   Poisson's    ratio   of    the   supporting  medium   is   changed   from  0   to  y, 
and   that    the   critical   buckling   pressure   for   frictionless    interface 
conditions    is   only   about    70%  as   large  as    the   critical   pressure   for   fully- 
bonded   interface   conditions.      Thus,    Forrestal   and   Herrmann's    theory 
indicates   that    slip   at    the   conduit-soil    interface   decreases   the  critical 
buckling  pressure   significantly. 

Luscher    (1966) 

In   his   study   of   buckling  of   soil-surrounded   tubes,    Luscher    (1966) 
formulated  a   semi-empirical    solution    for   the   critical,    uniform,    radially 
applied   buckling  pressure  acting   on   the  wall   of   a   buried  conduit. 
Luscher 's   theoretical   analysis   considered  buckling  of  an   elastically 
supported   ring,    which  he  believed  modeled   the   behavior   of   a   long    tube 
with   radial   elastic    support.      The   elastic    support    in    Luscher' s   model 
was    identical   to   that    in   Cheney's    (1963)    representation,    and   consequently, 
the   equation   developed  by   Luscher  was    identical   to   eqn    (7). 

_    EI(n   -1)  s , 

Pc   -  3  +  -J-  ^^°> 

R  n   -1 

-2 
K     =   modulus   of    soil    support   =  k      R    (FL      ) 

Equation    (10)   was   then   minimized  with  respect   to   n,    and   an   expression 
for   the   lower-bound  critical   buckling  pressure  was   obtained. 
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Pc=    ^ 


K      EI/R' 
s 


(lOi) 


Luscher's    treatment    varied   from  Cheney's    in   that   an   empirical   expression 
was   proposed   for   the  modulus   of    soil    support   based  on    the   results   of 
buckling   tests   on   1.63   inch   diameter   tubes   and  oedometer   tests   on   the 
supporting   sand.      The   relationship   between    the  modulus   of   soil    support, 
K   ,    obtained   from  buckling   tests,    and   the   constrained  modulus,    M   , 
obtained   from  oedometer   tests,    was   given   as 


B(l+v   )(l-2v   ) 

K     =    BE     =  r^ ^  M 

s  s  1-v  s 


(11) 


1   - 


R.^2 

_i  I 

R 


(1+v^) 


1+ 


R.^2 
o^ 


(12) 


R.  =  inside  radius  of  elastic  ring  of  soil  support 

R  =  outside  radius  of  elastic  ring  of  soil  support 

E  =  Young's  modulus  of  soil 
s       * 

M  =  constrained  soil  modulus 
s 


Comparison  of  results  from  tests  conducted  by  Luscher  to  predictions 
of  critical  buckling  pressure  by  eqn  (8)  showed  that   if  the  'experimental' 
values  of  n  are  used,  eqn  (8)  predicted  critical  pressures  much  larger 
than  those  measured  while  eqn  (lOi)  accurately  predicted  the  buckling 
loads.   Agreement  of  the  experimental  results  with  eqn  (lOi)  is  not 
surprising,  because  the  expressions  for  K  were  derived  from  the 
experimental  results.   Luscher  recognized  this  and  mentioned  that  eqn 
(lOi)  may  not  be  accurate  for  some  ranges  of  tube  and  soil  properties 
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because;  (1)  inelastic  buckling  was  not  considered,  (2)  experimental 
conditions  may  have  placed  artificial  constraints  on  the  tube's 
behavior,  and  (3)  the  test  results  may  not  be  scalable  to  larger  tubes. 

Based  on  his  studies,  Luscher  believed  that  elastic  buckling  rather 
than  wall  crushing  might  be  the  controlling  mode  of  failure  for  thin- 
walled,  smooth  metal  tubes  even  when  supported  by  dense  soil.   However, 
failure  would  occur  in  a  higher  buckling  mode  and  at  a  pressure  much 
larger  than  that  predicted  by  buckling  models  without  circumferential 
soil  support. 

Duns  and  Butterfield  (1971) 

In  developing  theoretical  stability  equations  for  long  cylindrical 
shells  in  an  elastic  medium,  Duns  and  Butterfield  used  many  of  uhe  same 
assumptions  employed  by  Forrestal  and  Herrmann  (1965)  except  that: 

(1)  nonlinear  terms  in  the  stability  equations  were  neglected 

(2)  displacements  associated  with  initial  stresses  were  not  considered 

(3)  buckling  in  high  modes  was  required  (n  >_  7). 

Essentially,  Bodner's  solution  was  extended  to  include  the  effects  of 

elastic  support.   As  in  previous  theories,  a  uniform  constant-direction 

pressure  bearing  on  the  conduit  produced  a  state  of  pure  hoop  compression 

in  the  conduit  wall.   When  the  thrust,  F  ,  reached  a  critical  value, 

c  ' 

given  by  eqn  (13),  failure  would  occur  by  elastic  buckling  in  high  modes. 


F  =   ^^ 


(1-V  ) 


+  k 
s 


2 

(13) 


F     =   critical   compressive   thrust    (FL      ) 
kg   =   modulus   of   soil   reaction    (FL      ) 
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Duns  and  Butterfield  took  advantage  of  their  simplifying  assumptions  to 
derive  a  theoretical  expression  for  the  modulus  of  soil  reaction,  k  ,  as  a 
function  of  the  soil  properties,  the  radius  of  the  conduit,  and  the  mode 
of  buckling,  given  by 

E^(n^-l) 

\^   R(l+V   ){2n+l-2v    (n+D)      ^^°''^^'^  soil-pipe    interface)  (14a) 

s  s 


Equation    (14a)    was   determined   by   means   of   the  Michell-Airy   stress   and 

displacement    functions.      An   approximate   expression    for   the  modulus   of 

soil    reaction   at   a    frictionless   soil-pipe   interface   was   also   derived   by 

minimizing   the  potential   energy,    giving 

k 
k'    =  — T-      (frictionless   interface)  (14b) 

The  mode  corresponding   to    the  critical   buckling  pressure   is   determined 
using   the  appropriate  eqn    (14)    and   eqn    (15)    below, 

2             2        ^s^^-^^)    ^^ 
^    ^^+1^    =— li (15) 

The  k^  modulus  is  initially  determined  (from  eqn  14)  for  an  assumed  n  and 
then  substituted  into  eqn  (15)  to  calculate  the  buckling  mode.   The 
iterative  procedure  ends  when  the  assumed  value  equals  the  buckling  mode 
number  calculated  from  eqn  (15). 

Duns  and  Butterfield  also  obtained  an  equation  relating  the  critical 
compressive  thrust  in  the  conduit  to  a  uniform  pressure,  p  ,  applied 
over  the  soil  surface.   This  was  accomplished  by  using  an  empirical 
expression  derived  from  finite  element  analyses  that  related  the 
compressive  thrust  to  the  surface  soil  pressure. 
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F  = 
c 


3(1-V  )  R  p 

S        V 

27v 


(16) 


This  expression  does  not  consider  the  conduit  stirfncss  and  appears  to 
be  valid  only  for  the  case  of  neutral  arching.   It  differs  from  eqa  (3) 
(White  and  Layer,  1960)  in  that  it  includes  an  effect  due  to  Poisson's 
ratio  of  the  soil. 

A  depth  of  burial  influence  factor,  similar  to  the  B  factor 
developed  by  Luscher  (eqn  12),  was  also  incorporated  into  the  analysis. 


1  - 


■R.\2 

1 

R 


\  = 


R.^2 

'•  0-' 


1  + 


(16a) 


F,    =   burial   depth   factor 

Substituting  eqns    (14a)    and    (16)    into   eqn    (13),    assuming  n    >   >  1,    and 
minimizing  with   respect   to  n,    the   lower-bound   critical    soil   over- 
pressure   is    given   by 


2-V          F^   '    E 
p      =  -ji — --^; 5—   (bonded   interface) 


V        (1-v   )    ,     ,,      2, 
s      4n(l-v   ) 


(17) 


where   the  critical   buckling  mode   is    given   by, 


F^E^R^l-V^) 

4(l-v   )EI 
s 


1/3 


(17a) 
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While  their  analysis  assumed  a  constant  compressive  thrust  in  the 
conduit  wall,  Duns  and  Butterfield  believed  that  there  was  little 
difference  between  the  buckling  stress  of  a  conduit  under  uniform 
load  and  one  under  non-uniform  load,  as  long  as  the  maximum  thrust 
developed  in  the  second  case  was  equal  to  the  uniform  thrust  in  the 
first  case.   The  development  of  bending  stresses  under  non-uniform 
loads  (or  of  initial  deflections  to  relieve  these  stresses),  which 
could  significantly  reduce  the  critical  buckling  stress,  was  overlooked. 

Chelapati  and  All  good  (1972) 

In  order  to  represent  the  boundary  conditions  at  a  culvert-soil 
interface  more  closely,  Chelapati  and  Allgood  (1972)  attempted  to  derive 
a  more  rigorous  theoretical  solution  for  buckling  of  deeply  buried 
culverts  under  a  uniform  external  pressure,  and  to  verify  the  solution 
with  extensively  instrumented  model  tests.   The  theoretical  solution, 
using  a  bonded  interface  and  minimum  total  potential  energy  expressions 
from  large  deflection  theory,  treated  two  cases  of  soil  support.   The 
first  case  assumed  continuous  elastic  radial  support  around  the  entire 
circumference,  and  the  second  case  involved  soil  support  only  on  outward 
deflecting  lobes  of  circumferential  buckling  waves,  thus  decreasing  the 
critical  buckling  pressure.   Luscher's  (1966)  expression  for  K   {eqn 
(11)}  was  used  directly  by  Chelapati  and  Allgood  to  describe  the  soil 
support  properties;  however,  Chelapati  and  Allgood  applied  the  analysis 
to  a  shell  rather  than  a  ring.   The  resulting  general  solution  for  the 
critical  uniform  external  pressure  acting  on  a  culvert  with  continuous 
circumferential  support  was  given  as 
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(i^^im  ^  _  J^  (^3, 


'  c       2   3     2 

2 
which  is  identical  to  eqn  (10),  except  for  the  (l-v'")  term.   In  the  case 

where  soil  support  is  only  on  outward  deflecting  lobes,  K  in  eqn  (18) 

is  halved.   Minimizing  eqn  (18)  with  respect  to  n  and  assuming  V  =  0.3 

gives 


P.  =  2.1 


] 


K  EI/R^  (18i) 


The  corresponding  expression  for  the  lower  bound  buckling  pressure  of  a 
culvert  supported  by  soil  on  outward  deflecting  lobes  only  is  given  by 


P,  =  1.5 


\. 


EI/R  (18ii) 


It  is  to  be  noted  that  eqns  (18i)  and  (ISii)  are  lower  bound  solutions 
and  should  be  expected  generally  to  yield  conservative  values  for  the 
critical  buckling  pressure.   Indeed,  eqn  (18i)  was  found  to  give  a 
conservative  prediction  of  the  results  from  small  scale  laboratory  tests 
conducted  by  Allgood  and  Ciani  (1968). 

Cheney  (1976) 

Cheney    (1976)    revised   his   derivation  of   eqn    (7)    and   obtained   an 
expression    for  buckling  of   flexible   tubes   under  a   uniform  external  pressure 
and   supported  by  an   elastic   continuum.      He  derived  an   expression    for 
k     which  previously   had   been  obtained  approximately    (Duns  and   Butterfield) 
or   determined   semi-empirically    (Luscher).      The  resulting   form  of   this 
expression   is   given  by, 
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k     = 
s 


E^(n   -1) 


R(l-v   )(2n+l-v  ) 
s  s 


(19) 


E     =    secant   elastic   modulus   of   soil 
s 


Cheney's   expression    for   the   critical   uniform  buckling  pressure   is    given   by 


EI    (n   -1)   ^ 
p^  =  ^  '^  + 


R^   (1-v   )         (2n+l-v   )(l-v^) 
s  s 


(20) 


which  is  then  minimized  with  respect  to  n  to  give. 


P,  =  1.2 


1/3 

EI 

2      3 
(1-V'')R^ 

(1-v^) 
s 


2/3 


(20i) 


where  p^  in  eqn  (201)  is  the  lower  bound,  critical  buckling  pressure.   However, 
when  this  equation  was  used  to  predict  the  test  results  obtained  by  Allgood 
and  Ciani  (1968),  it  made  unconservative  predictions.   Cheney  believed 
that  soil  representation  by  the  elastic  secant  modulus  caused  this 
overprediction  of  critical  pressure,  and  he  proposed  the  use  of  a  La  igent 
elastic  soil  modulus  related  to  the  stress  state  in  the  soil  arouiiL.  Che 
tube  to  provide  a  more  valid  predictive  model.   Cheney  suggested  the 
use  of  'Critical  State  Soil  Mechanics'  to  estimate  a  tangent  modulus 
value,  but  these  concepts  were  not  sufficiently  developed  for  practical 
use  at  the  time. 

By  reducing  the  modulus  E   in  eqn  (20i)  to  1/8  the  secant  value 
eqn  (21)  resulted,  vjhich  fit  the  experimental  results  of  Allgood  and 
Ciani  quite  well.   Cheney  concluded  that  the  effective  tangent  modulus 
must  be  one-eighth  the  value  of  the  secant  modulus   and  proposed 
eqn  (21)  as  a  design  formula. 


23. 


p  =  0.30 
c 


EI 


3    2 


1/3 


s 


2/3 


(21) 


The  validity  of  Cheney's  empirical  adjustment  of  E  remains  to  be 

proven  by  further  tests  in  which  the  soil  stress  state  at  the  onset  of 

buckling  is  known  and  a  corresponding  tangent  modulus  can  be  determined 
independently . 

Meyerhof  and  Baikie  (1963) 

Meyerhof  and  Baikie  extended  plate  buckling  theory  developed  by 
Timoshenko  and  Gere  (1961)  to  approximate  the  buckling  behavior  of 
elastically  supported  cylindrical  shells.   This  approach,  although 
developed  before  the  more  rigorous  shell  buckling  theories  that 
considered  soil  support  were  formulated,  is  used  in  current  design 
practice. 

In  the  theoretical  development,  solutions  for  buckling  of 
long  flat  plates  supported  by  a  Winkler  foundation  (independent  springs) 
with  a  constant  subgrade  modulus  were  combined  with  the  assumption 
that  the  critical  buckling  shape  of  a  cylindrical  shell  approaches 
that  of  a  flat  plate  at  very  high  buckling  modes.   Thus,  eqn  (2) 
developed  by  Timoshenko  and  Gere  was  combined  with  the  equation  for  the 
critical  buckling  force,  F  ,  of  an  elastically  supported  long  flat 
plate  to  obtain. 


F  = 
c 


EI 


2   2 
(l-V^)r 


2    4  ' 
(1-v  )k  R 

(n+1)^  -  1  + ■-— 

{(n+1)  -1}EI 


(22) 


k  =  subgrade  modulus 
s      ° 
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The  suberade  modulus,  k  ,  was  then  related  to  the  soil  modulus,  E  ,  and 

s  s 

Poisson's   ratio,    V   ,    by    the   following  approximate    relationship, 


k     = 
s 


2(l-v   )R 
s 


(23) 


Substitution  of  eqn  (23)  into  eqn  (22)  yields. 


F  = 
c 


EI 


(1-v^) 


f(n+l)  -1)1 


(1-v^)  2{(n+l)^-l} 
s 


(24) 


which  reduces  approximately  to  eqn  (13)  developed  by  Duns  and  Butterfield 
for  large  values  of  n. 


Equation  (24)  was  minimized  with  respect  to  n  and  divided  by  the 
cross-sectional  wall  area  to  obtain  an  expression  for  the  critical 
buckling  stress,  f  ,  valid  for 


3„    15 


R  E 


(25) 


and  given  by. 


^c4 


EI  k 
s_ 

(1-v^) 


(26) 


where,  A  =  x-sectional  area  of  the  wall 
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Eqn  (26)  is  equivalent  to  eqn  (181)  developed  by  Chelapati  and 
Allgood. 

An  empirical  modification  of  eqn  (26)  was  then  made,  similar  to 
the  Rankine-Gordon  equation  developed  for  column  buckling,  to  extend 
its  applicability  to  inelastic  buckling  of  less  flexible  conduits. 
The  resulting  expression. 


f 


(l-V^)(l-vJ)R 


(27) 


2  E   EI 
s 


where  f  =  yield  stress 

y 

gives  much  lower  critical  stresses  than  eqn  (26),  and  it  was  considered 
that  the  additional  safety  factor  could  compensate  for  possible 
conduit  imperfections.   However,  the  validity  of  the  empirical 
adjustments  has  not  been  verified  for  conduits  with  soil  support. 

Meyerhof  (1968)  suggested  adjustment  of  eqn  (2  7)  for  analyzing 
culverts  with  depth  of  cover  less  than  one  diameter  by  using  a  reduced 
soil  modulus  given  by  Luscher's  eqn  (11)  and  multiplying  the  critical 
buckling  stress  by  a  reduction  factor  n  accounting  for  the  loss  of 
soil  support  due  to  shallow  cover,  where 


n  =  4  (28) 


0.1<2|<1.0 


where  z  =  depth  of  soil  cover  above  crown 
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B.   Elastic  Local  Buckling 


Amstutz  (1953) 


Amstutz  (1953)  presented  a  solution  for  the  critical  uniform 
external  hydrostatic  pressure  producing  a  local  inward  buckle  in  a 
flexible  cylindrical  shell  rigidly  restrained  from  outward  buckling. 
Amstutz  considered  bending  stresses  as  well  as  thrust  stresses  that 
developed  in  the  shell  wall  during  buckling.   Calculation  of  the 
bending  stress  along  the  buckled  section  was  accomodated  by  defining 
the  buckling  displacements  as  a  cosine  function.   From  the  known 
buckled  shape  the  radius  of  curvature  and  the  corresponding  bending 
stress  could  be  obtained  at  each  point.   The  extreme  fiber  stress  was 
not  allowed  to  exceed  the  yield  stress  of  the  shell  material  and, 
thus,  the  local  buckling  criterion  was  defined.   The  critical  normal 
stress  was  determined  by  the  following  lower  bound  (n  =  1)  solution. 


f^(l-v-) 


1  +  12 


2  f  (1-v^) 
c 


_  2R  (V'c^^l-^  ) 


L 


t        E 


(29 


where  f  =  yield  stress 

y 

f  =  critical  circumferential  thrust  stress 
c 

This  critical  thrust  stress  can  be  related  to  the  critical  external 
hydrostatic  pressure  by  eqn  (3). 


Cheney  (19  71) 

Cheney  (1971)  reviewed  Luscher's  analysis  and  discussed  some 
unsatisfactory  aspects  of  his  results.   Luscher's  K  values  did  not 
seem  to  reflect  effects  of  soil  density  on  the  supporting  capacity  of 
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the   soil.      Back-calculation  of   K     for  loose  and   dense   sand   from  the 

s 

model  buckling  tests  should  have  yielded  different  values  but,  in 
fact,  they  were  found  to  be  essentially  the  same.   Cheney  believed 
this  discrepancy  was  probably  due  to  the  in-placc  density  of  Llie 
sand  in  the  tests  changing  to  a  value  between  the  dense  and  loose 
states  before  the  full  load  was  applied  to  the  soil,  and  thus  the 
soil  density  at  buckling  was  essentially  the  same  for  all  the 
laboratory  tests.   However,  such  behavior  may  not  have  been  possible 
in  the  small  scale  at  which  the  tubes  were  tested,  because  the 
strains  associated  with  buckling  of  small  tubes  are  too  small  to 
induce  appreciable  changes  in  sand  density. 

Cheney  also  pointed  out  some  conditions  in  the  tests  conducted  by 
Luscher  that  may  not  have  been  properly  modeled  by  his  theory. 

(1)  External  pressures  exerted  on  the  tubes  originated  from 
two  sources;  pore  pressures  from  water-filled  voids  of 

the  soil  cylinder  surrounding  the  tube  and  effective  stresses 
from  the  soil.   These  two  types  of  load  represented 
hydrostatic  and  non-hydrostatic  pressures,  respectively, 
while  Luscher' s  theory  only  considered  constant-direciton 
loads . 

(2)  The  tubes  tested  by  Luscher  generally  failed  after  the  formation 
of  a  snap-through,  local  buckle.   Theoretically,  failure  was 
assumed  to  take  place  in  high  buckling  modes  involving  the 
entire  tube  circumference.   Thus,  Luscher' s  theory  was  not 
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consistent  with  the  observed  buckling  mode,  which  strongly 
mitigates  against  its  general  applicability.   Also,  the  effect 
of  local  inward  displacements  on  the  pressure  distribution 
around  the  tube  was  not  considered  in  Luscher's  analysis. 

Cheney  attempted  to  model  theoretically  the  local  buckling  behavior  of 
a  buried  conduit  subjected  to  an  initial  uniform  radial  pressure.   As 
illustrated  in  Figure  3,  he  assumed  rigid  boundary  restraints  preventing 
outward  buckling,  as  was  the  case  in  Amstutz's  theory.   However,  a  solution 
for  the  critical  uniform  external  radial  pressure  initiating  a  loc.il  buckle 
was  obtained  by  minimizing  the  energy  of  distortion  into  the  buckled  shape 
to  determine  the  governing  equilibrium  condition.  The  total  energy  of 
distortion  was  given  as  the  sum  of  the  internal  strain  energy  due  to 
extension  and  bending  minus  work  done  by  the  boundary  forces  and  moments. 
Unlike  Amstutz's  model,  in  which  the  boundary  hydrostatic  pressure  remained 
constant  during  buckling,   Cheney  attempted  to  represent  the  true  behavior 
of  the  soil  at  the  location  of  inward  deflection.   As  the  deflection 
increased  the  boundary  pressure  decreased  according  to  a  linear  force- 
deflection  relationship  (Figure  3),   A  graphical  procedure  was  developed 
to  obtain  solutions  to  the  complex  differential  equations  applicable  to 
these  special  boundary  conditions.   The  final  shape  of  the  symmetric  buckle 
was  determined  by  the  equilibrium  equations  and  was  not  arbitrarily  taken 
to  be  a  specific  cosine  function  (as  assumed  by  Amstutz). 

Cheney  noted  that  Luscher's  soil  support  modulus  varied  experimentally 
with  the  thickness  of  the  supporting  soil  ring  and  the  effective  soil 
pressure  producing  buckling.   He  subsequently  developed  an  empirical 
expression  for  K  given  as. 
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FIGURE  3         DISTRIBUTION  OF  CRITICAL  PRESSURE   DURING 
LOCAL  BUCKLING  OF  SOIL- SURROUNDED  TUBES 
(after  CHENEY.  1971) 
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K     =    720   a' 
s 


R   -R. 
o       1 

R. 

1 


0.153 


(30) 


where,    a'  =   effective   soil   pressure   at    the   conduit    soil    interface 

R    ,R.    =   outside   and   inside    radii  of    the    ring  of   soil   support 

O  1  b  Kf 

surrounding   the   conduit 
When   the   results    from  Cheney's   analysis   were   compared  with   Luscher's 
experimental   data   they   overestimated   the  observed  critical   buckling 
pressures.      Although  a   rigid   boundary   condition   appears   severe,    Cheney 
blamed   the   disagreement   on   imperfections    in    the  loaded   tubes   that   were 
not   accounted   for   theoretically   and  noted   that   the   trends   of   buckling 
behavior  predicted   by   his   theory   paralleled   the   trends   shown   by   the 
experimental    results.      By   using  an   effective   wall    thickness    to   account 
for   structural    imperfections,    Cheney   was   able   to    fit    his   theoretical 
results   to   the   experimental   values.    However,    Cheney   conceded   that   more 
research  was   needed   before   his    theory   could   be   used   practically   to   predict 
local   buckling.      He   recommended   full-scale   laboratory   testing  of   conduits 
to   obtain   more   realistic  notions   of    the   importance  of   local    structui^al 
imperfections    in   the   pipe.      Also,    results    from  full-scale   bucklin,^   tests 
could  verify   the   validity   of   the   soil    support   modulus   expressions 
obtained   from  small-scale   tests. 


El-Bayoumy    (1972) 

El-Bayoumy   analyzed   the   critical   uniform  pressure  producing   snap- 
through   buckling   in   a   rigidly   confined,    flexible   circular   ring  with 
boundary   conditions    identical    to    those   in   the  problem  treated   by   Cheney 
(1971).      He  also   used  energy   methods   to   obtain   the   equilibrium  equations 
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for  the  unbuckled  configuration.   However,  the  boundary  pressure  was  not 
reduced  at  the  location  of  buckling.   Rather,  a  snap-through  buckle  was 
assumed  to  occur  when  the  critical  pressure  was  reached,  and  the  rapidity  of 
the  mechanism  was  such  that  a  reduced  pressure  at  the  buckle  location  did 
not  affect  the  buckling  energy.   Differential  equations  were  developed 
using  variational  calculus  which  related  the  strains  and  rotations  of  the 
buckled  configuration  to  the  circumferential  stress  in  the  ring. 
Solution  of  these  equations  gave  lower  bound  buckling  loads  identical  to 
those  obtained  by  Cheney  (1971)  for  the  constant  pressure  case  (k^  =  0). 
However,  the  detailed  development  revealed  some  important  energy  considerations 
clarifying  the  snap-through  buckling  phenomenon. 

As  an  increasing  external  pressure  is  applied  to  a  rigidly  confined 
ring,  the  ring  contracts  in  hoop  compression  and  the  stored  strain  energy 
increases.   When  the  strain  energy  reaches  the  level  corresponding  to  the 
final  buckling  state  the  ring  has  a  tendency  to  buckle  and  reach  a  lower 
energy  state,  but  in  order  to  do  so  it  must  pass  through  an  initial,  small 
deflection  buckling  phase  which  exists  at  a  higher  energy  level.   An 
external  disturbance  or  a  ring  imperfection  is  required  to  overcome  the 
energy  barrier  that  otherwise  prevents  the  initiation  of  buckling.   If  such 
a  disturbance  is  provided,  the  ring  passes  rapidly  from  the  unbuckled 
to  the  initial  buckling  state,  and  then  to  large  deflection  buckling,  with 
an  accompanying  sudden  release  of  energy.   To  obtain  his  lower-bound  solution 
for  the  critical  external  pressure,  El-Bayoumy  assumed  that  an  imperfection 
existed  which  would  allow  snap-through  buckling  to  occur  when  the  strain 
energy  of  the  unbuckled  ring  equalled  that  of  large  deflection  buckling 
equilibrium.   Theoretically,  snap-through  buckling  will  not  occur  in  perfect 
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rings  free  from  external  disturbance.   Consequently,  flexible  conduits 
without  imperfection  subjected  to  uniform  external  pressure  should  ♦'ail 
by  some  other  mechanism  such  as  elastic  buckling  in  high  modes.   The 
fact  that  snap-through  buckling  occurs  more  frequently  in  practice 
indicates  that  geometric  and  material  imperfections  (and  possibly 
residual  stresses)  play  an  important  role  in  conduit  performance. 

Kloppel  and  Clock  (1970) 

Extensive  investigations  of  buried  flexible  conduits  by  Spangler 

(1941)  led  to  the  conclusion  that  vertical  deflections  on  the  order  of 

20%  of  the  diameter  would  generally  induce  an  instability  in  a 

conduit.   Spangler  reasoned  that  if  the  deflections  were  limited  to 

5  percent  a  substantial  safety  factor  against  instability  would  be 

provided,  and  this  criterion  has  been  used  extensively  in  culvert  design, 

Kloppel  and  Clock  (1970)  set  out  to  develop  a  rational  method  of 

predicting  the  instability  of  initially  deflected,  buried  conduits. 

They  distinguished  two  zones  of  interaction:  one  located  in  the  upper 

portion  where  the  conduit  wall  deflects  away  from  the  surrounding  soil 

and  is  subjected  to  "active"  earth  pressure,  and  the  other  on  the 

lower  portion  where  the  conduit  wall  presses  into  the  surrounding  soil 

and  is  subjected  to  "passive"  earth  pressure.   Because  of  the  reduced 

soil  support  of  the  upper  section,  it  was  reasoned  that  an  instability 

is  most  likely  to  occur  near  the  crown  of  the  conduit.   Kloppel  and 

Clock  chose  to  model  the  upper  conduit  section  by  a  hinged  arch  with 

"     E    "^ 


radial  elastic  (soil)  support 


R(l  +  V  ) 


along  its  circumference  and 


tangential  and  rotational  elastic  support  (k  and  k^)  at  the  hinges, 
thereby  representing  soil  support  and  restraint  provided  by  the  lower 
conduit  section.   The  ensuing  analysis  considered  the  following  conditions: 
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(a)  Eliiptical   as   well   as   circular   conduits 

(b)  Friction   at   the   conduit-soil    interface 

(c)  Nonuniform   radial   pressures   having   a   maximum  value   at    the 
crown    (so    that   buckling  was   always    initiated  at    the   crown) 

(d)  Modulus   of   soil   restraint   either   constant   or   stress-level 
dependent 

(e)  Symmetric   displacements   prior   to   buckling 

(f)  The   influence  of   a   plastic    hinge  at    the   crown. 
Nonlinear  terms    in    the   stability   expressions   were   retained  and 

energy   methods   employing   variational   calculus   were   used   to   analyze   the 
following   three   representations   of   potentially   unstable   flexible 
conduits: 

(1)  two-hinged  arch  with   rotational  and  tangential   restraint  at 
the   hinges, 

(2)  two-hinged  arch  neglecting   rotational    restraint   at    the  hinges,    and 

(3)  three-hinged  arch    (one  hinge  at   the  crown)   with  no   rotational 
restraint   at   the   hinges    (to   model   the   buckling   resistance   of 
a  conduit   after   formation  of   plastic   hinges). 

In   each  analysis   a  nonuniform   distribution   of   radial    "active"   soil 
pressure  was   expressed  as   a   function   of   the   coefficient   of   active 
earth  pressure,    A.      The  horizontal   active   pressure   at   the   springline   was 
defined  by, 


P         .        =   p      A  (31) 

spring  s 


where  p   is  the  vertical  crown  pressure  equal  to  the  sum  of  the  overburden 


and  live  load  pressures  at  the  crown. 
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The  distribution  of  active  radial  pressure  about  the  conduit  circum- 
ference was  expressed  as  a  cosine  function  (Fig.  4A) 


p  =  p   cos 
s 


(32a) 


^  = 


4  cos   (X) 


(32b) 


In    formulating   the   distribution,    Kloppel   and   Clock,  neglected  pressure 

redistribution   due   to   soil   arching  around   the   conduit.      A  value  of   0.5 

3tt 
was   commonly  assigned   to    A,    which   set    i|j     equal   to  —j-  radians. 

In  order   to    solve   the   instability   problem,    an   arch   section  of   the 

conduit    (defined   by   2    c{)      in   Fig.    4B)   was   analyzed.      Restraint   provided 

by   the   remainder  of    the  conduit   was   represented  by    rotational   and 

tangential   elastic   moduli,    k     and   k    ,    at    the   arch   supports.      The 

R       1 

boundary  pressure  distribution  on  the  arch  section  was  decomposed  into 
a  uniform  component  p  ,  and  a  nonuniform  component  p, ,  so  that 


P  =  Pq  +  Pj^  sin 


(33) 


where,  p  =  p   cos  -z — — 
^o   ^s     2  (jj 


and  p  =  p  +  p  =  vertical  overburden  plus  live  load  pressures  at 


To  solve  the  snap-through  instability  problem  an  in  It  Lai  value  of 
A  was  chosen  and  deformation  into  a  single  mode  buckle,  defined  by  a 
cosine  function,  was  considered.   A  load-deformation  expression  relating 
the  thrust,  F,  in  the  conduit  wall  at  the  crown  to  the  amplitude  of  the 
deflection  was  developed  by  minimizing  potential  energy.   A  second 
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FIGURE  4A   ACTIVE  PRESSURE  DISTRIBUTION  AROUND 
CIRCULAR  CONDUIT 


FIGURE  4B  PRESSURE  DISTRIBUTION  ON  TWO-HINGED 
ARCH  (after  KLQpPEL  AND  CLOCK,  1970) 
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expression  relating  F  to  the  crown  pressure  p   for  a  given  boundary 
pressure  distribution  was  obtained  from  equilibrium  requirements. 
Values  of  F  were  then  assumed  and  the  corresponding  deflected  shape 
and  p   were  calculated.   The  critical  thrust  F  for  the  chosen  (f)  was 

^S  CO 

determined  when   an   increment    in   F  produced  no    increment    in   p    .      At 

^s 

this  point  increasing  the  load  did  not  increase  the  potential  energy 

of  the  conduit-soil  system,  and  an  instability  was  considered  to  occur. 

Other  values  of  d)  were  chosen  and  the  minimum  thrust  at  which  an 
o 

instability  occurred  was  taken  to  be  the  critical  buckling  load. 

A  similar  analysis  was  performed  for  high  mode  buckling,  except 
that  instead  of  restricting  deformations  to  a  single  snap-through 
buckle,  deflection  into  a  large  number  of  waves  defined  by  a  periodic 
trigonometric  function  was  allowed.   The  subsequent  expressions  for  the 
critical  buckling  load  were  minimized  with  respect  to  the  buckling  mode 
and  a  high  mode  buckling  expression  similar  to  the  one  derived  by  Luscher 
(eqn  lOi)  for  a  flexible  ring  was  obtained. 

Sensitivity  analyses  were  performed  by  Kloppel  and  Clock  to  determine 
the  influence  of  the  assumed  boundary  conditions  on  the  calculated  snap- 
through  buckling  load,  which  led  to  the  following  general  conclusions 
(see  Fig.  5)  : 

(1)   For  small  values  of  the  conduit-soil  flexibility  factor, 
EI(l+v^)       _^ 

F  =  — ~ <  10   ,  variation  of  the  conduit  or  soil 

^    r'  E 
stiffness  nas  little  effect  on  the  snap-through  buckling 

load. 
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FIGURE  5      VARIATION  OF  CRITICAL   THRUST  FORCE   WITH 
BOUNDARY  CONDITIONS   AND  BUCKLING  MODE 
(after   KLOPPEL  AND   GLOCK,  1970) 
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(2)  Neglect  of  rotational  restraint  of  the  arch  support  (k  ) 

R 

reduces   the  critical   buckling   load. 

(3)  Given    identical   conduit   and   soil   stiffness   and   equivalent 
boundary   pressure   distributions,    the  critical   buckling   load 
of   an   elliptical   conduit,    with   the  major  axis   horizontal,    is 
approximately   one-half   that   of   a  circular  conduit. 

(4)  The  amount   of   soil-conduit    interface    friction   has  negligible 
influence   on    the  magnitude   of   the   critical    buckling   load. 

(5)  Increasing   the  nonuniformity   of   the  active   pressure   distribution 
(Figure  4A)    decreases   the   critical   buckling   load. 

(6)  Consideration  of   a   stress-level    dependent   elastic    soil   modulus, 

which   is   constrained   to    increase   linearly  with   stress,    slightly 

increases   the  critical   load.       {in    fact,    E     may   decrease  with 

increasing   deformation.      Kloppel   and   Clock  concluded   that 

because   the   stress-level    dependence  was   uncertain,    a   constant 

value  of    E      should   be  assumed. } 
s 

Comparison   of   the  critical   loads    in   high  mode   and   snap-through 

-4 
buckling   showed   that    for   F      -   10        the  critical   loads   are  nearly 

identical.      This    indicates   that   very    flexible  conduits   with   good   soil 

-4 
support   will   tend   to   buckle    in   high  modes.      However,    for   F      >  10 

the   critical   snap-through   buckling   load   is   much  lower   than    the  high 

mode   buckling   load,    so    that    in    this   range   high  mode   buckling   theory    is 

inapplicable    (see   Fig.    5).       Investigation  of   conduit    instability   after 

formation  of  a  plastic   hinge  at    the  crown   resulted   in   a   lower 

prediction  of   the  critical   load. 
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Although   the   theoretical   studies    conducted  by   Klopj^el    and   r.lock 
were   of  a  more   general  nature   than  previous  work,    there  still   remain 
oversimplifications   that  prevent  close  modeling  of  the  true  soil-conduit 
system.      Their   formulation  of  the  boundary  pressiore   distribution  mis- 
represents  the   true  distribution  around  buried   flexible   conduits  ,    and 
is   not   compatible  with   the   assumed  deflection  pattern  prior  to  buckling. 
Experience  has   sho-.^  that  the  normal  pressure   at   the   crown   is   often 
smaller  than   the  norainal  overburden   pressure,    and  the   springline 
pressure   is   larger  than   the   active  pressure,  because  of  arching  that 
accompanies    conduit   deflection.      Thus-,   the  actual  pressure   distribution 
is    less    conducive  to  buckling  at   the   cro-^n   than  the  assumed  one,    and 
Kloppel  and  Clock's    lover  bound  theory   may  be  a  reasonable  solution    ^o 
the  buckling  problem. 

C.      Inelastic   Buckling  Theory 
Inelastic  buckling  of  cylindrical   shells    is    a   complex  problem  that 
has   been    analyzed  only   for   simple  boundary   conditions    and   load   distributions. 
Several   investigators    (Bijlaard,   19^9;    Gerard,    1957;   Reynolds,    19^0;    Lee, 
1962)   have   examined  inelastic  buckling  of  cylinders    of  finite  length 
subjected  to   uniform  external  pressure  but  without   circumferential    '.Al 
support.      Gerard   (1957)    proposed   the   use   of  the    following  expression    for 
an   "effective"   modulus    to  be   substituted   into   expressions    derived   for 
elastic  buckling  to  predict  the   critical   inelastic  buck"'ing  pressure: 
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^eff 


E  E. 


1  -  V 


1  -  V 


(3A) 


where  E  =  secant  modulus 
E  =  tangent  modulus 

V  =  elastic  Poisson's  ratio 
e 

V  =  plastic  Poisson's  ratio 
P 

The  above  expression  is  not  applicable  if  E   =  0  after  yielding.   In  such 
cases  an  empirical  value  of  E  ^^  must  be  obtained.   In  addition,  the  use 
of  an  effective  modulus  may  be  applicable  only  to  materials  that  are 
free  of  residual  stresses,  such  as  heat-treated  aluminum  alloys.   If  an 
accurate  representation  of  the  elastic  and  plastic  stress-strain  behavior 
is  available  a  trial  and  error  approach  to  determination  of  the  effective 
modulus  may  provide  a  good  estimate  of  the  critical  inelastic  buckling 
stress — as  long  as  geometric  imperfections  are  small.   However,  in  cases 
where  residual  stresses  are  significant  the  use  of  an  effective  modulus 
is  inapplicable. 

A  simple,  yet  reasonable,  approach  to  prediction  of  critical  loads 
in  the  inelastic  range  was  proposed  by  Gjelsvik  and  Bodner  (1962). 
Assuming  elastic-purely  plastic  behavior,  they  suggested  the  use  of 
elastic  instability  equations  on  the  unyielded  fraction  of  the  section 
to  predict  the  onset  of  buckling.   Thus,  an  equivalent  section  was 
proposed  rather  than  an  effective  modulus,  which  may  be  applicable  to 
less  flexible  conduits  subjected  to  relatively  large  bending  stresses. 

No  theories  have,  as  yet,  been  advanced  to  predict  inelastic 
buckling  of  radially  supported  conduits. 
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II.   COMPARISON  OF  ELASTIC  BUCKLING  THEORIES 

Many  formulas  used  for  predicting  critical  elastic  buckling  loads 
in  flexible  conduits  have  been  reviewed  in  the  previous  section.   They 
include  a  wide  variety  of  assumptions  and  boundary  conditions  that 
should  be  considered  when  evaluating  their  applicability  to  the  prediction 
of  buried  conduit  performance.   Table  1  highlights  the  significant 
assumptions  and  boundary  conditions  inherent  in  the  respective 
theoretical  formulations. 

The  review  of  theoretical  developments  in  Section  I  includes  a 
number  of  observations  and  comparisons.   For  convenience,  the  main 
points  are  summarized  below: 

a)  Analysis  of  shells  is  more  complicated  than  that  of  rings 
yet,  if  all  other  assumptions  are  comparable,  the  predicted 
critical  pressures  do  not  differ  greatly,  with  the  ring 
analysis  being  a  little  more  conservative. 

b)  Critical  pressures  that  remain  normal  to  the  conduit  wall 
are  slightly  lower  than  those  that  are  maintained  in  a 
constant  (radial)  direction. 

c)  Duns  and  Butterfield  related  the  critical  buckling  mode  to 
conduit  and  soil  properties.   When  these  critical  buckling 
modes  are  inserted  in  the  other  theories,  the  theoretical 
predictions  are  much  higher  than  the  lower  bound  values. 
This  is  because  all  theories  initially  assumed  buckling  in  a 
high  mode,  but  only  Duns  &  Butterfield  determined  the 
critical  high  mode  (for  n  ^  7).   Other  theories  are  minimized 
with  respect  to  the  mode  number  n,  regardless  of  whether  or 
not  the  resulting  mode  is  compatible  with  the  conduit  and 
soil  properties. 

d)  The  effect  of  slip  at  the  soil-conduit  interface  is  not 
definitely  known.  Two  theories  (Forrestal  and  Herrmann 
and  Duns  and  Butterfield)  show  that  slip  decreases 

the  critical  pressure,  but  they  disagree  considerably  on 
the  amount  of  reduction.   Kloppel  and  Clock,  on  the 
other  hand,  concluded  that  the  reduction  is  negligible. 
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(e)  An  important  source  of  disagreement  among  the  theories 
predicting  elastic  buckling  is  the  formulation  of  the 
soil  support  modulus.   Because  of  the  assumptions  and 
simplifications  characterizing  each  analysis,  several 
different  expressions  for  the  soil  support  modulus  have 
evolved.   Some  have  been  developed  by  minimizing  the 
critical  buckling  mode  so  that  a  lower-bound  solution  is 
obtained.   Others  have  been  similarly  derived  but  include 
empirical  factors  which  attempt  to  fit  the  theoretical 
predictions  to  observed  behavior.   The  soil  support 
moduli  in  various  theories  are  compared  in  Table  2, 

for  deep  soil  cover  and  a  soil  Poisson's  ratio  of  0.4. 
It  is  evident  from  this  comparison  that  lower-bound 
solutions  involve  large  reductions  in  the  soil  modulus 
not  necessarily  compatible  with  the  physical  aspects  of 
the  problem. 

(f)  With  the  exception  of  Kloppel  and  Clock  (1970),  all  theories 
deal  only  with  uniformly  applied  boundary  pressures  on 
circular  x-sections.   Nonuniform  pressures  adversely  affect 
the  critical  buckling  load,  and  elliptical  sections  buckle 
more  readily  than  circular  ones. 

(g)  Two  theories  (Donnell,  1956  and  Meyerhof  and  Baikie,  1963) 
consider  the  influence  of  geometric  imperfections  on  the 
critical  buckling  load;  both  theories  show  that  an  increase 
in  the  yield  strength  decreases  the  influence  of  a  given 
imperfection  level  on  the  critical  load. 

(h)   Four  theories  of  local  buckling  were  reviewed.   Three  of 
these  (Amstutz,  Cheney,  El-Bayoumy)  considered  the  snap- 
through  phenomenon  of  an  undeformed  conduit  rigidly 
restrained  from  outward  deflection  and  subjected  to  uniform 
pressure.   A  symmetric,  sinusoidal  snap-through  deflection 
was  forced  in  each  case.   Amstutz  considered  instability 
to  occur  when  the  extreme  fiber  stress  (including  bending 
stresses)  reached  the  yield  stress,  while  Cheney  and 
El-Bayoumy  arrived  at  the  instability  condition  through 
energy  considerations.   Partly  because  of  the  nature  of 
Amstutz 's  failure  assumption,  his  theory  is  a  lower  bound 
prediction.   With  no  reduction  in  pressure  during  buckling, 
results  from  Cheney  and  El-Bayoumy  are  identical.   Reduction 
of  pressure  as  the  conduit  snaps  through,  as  considered  by 
Cheney,  appears  to  be  unnecessary  because  of  the  nature  of 
the  energy  barrier  involved  in  the  snap-through  phenomenon. 
On  the  other  hand,  rigid  restraint  to  outward  deflection 
in  the  above  theories  tends  to  overestimate  the  buckling 
pressure. 

The  fourth  theory,  advanced  by  Kloppel  and  Clock,  considered 
nonuniform  boundary  pressures,  initial  deflections,  and 
variable  interface  friction.   The  results  showed  that  if  the 
soil-conduit  flexibility  factor  (F  )  is  low,  a  high  mode 


TABLE  2 

COMPARISON  OF  THEORETICAL  FORMULATIONS  OF 

THE  SOIL  SUPPORT  MODULUS 
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Investigator 


Suggested  Expression 
For  k 


Lower- Bound 

Value  of  k 
s 

(v^  =  0.4) 


Luscher  (eqn  11) 


r- 

1 

^■)^ 

(a) 

E 
s 

1     - 

^V 

rR.^2 
(1+V  ){l+hr^  (1-2V  )}R 


0.71  E 


Duns  &  Butterfield 
(eqn  14a) 


Cheney  (eqn  19) 


R(l+V  ){2n+l-2v  (n+1)} 
s         s 


R(l-V^)(2n+1-V  ) 
s        s 


4.0  E 


3.9  E 


(b) 


(b) 


Meyerhof  &  Baikie 
(eqn  23) 


2(1-V  )R 
s 


0.60  E 


Kloppel  and  Clock 
(1970) 


R(l+V   ) 


0.71   E 


(a)  f^il^ 

For  large   fill   heights,    R      >   >   R.    and     —       approaches    zero, 

*■  o^ 


(b) 


For  n   =    7. 
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buckle  was   more   critical    than   a   single  mode   snap-through, 
and  vice-versa    (Fig.    5).      Thus,    local    crimping   in    flexible 
conduits  with   good   soil   support   most   probably   originate   from 
significant    imperfections   or   local   residual    stresses    in    the 
conduit   wall.      Snap-through  buckling   is   also    influenced   by 
imperfections   and   local   residual    stress.      As   none  of    the 
buckling   theories   consider  all   the   influencing    factors   -   non- 
uniform pressure   distribution,    imperfections,    and  residual 
stresses   -   the   lowest   buckling  load   calculated   from  any   of   the 
revelant   theories   is   given  in  Fig.    6. 

The   lower-bound  critical   high  mode   buckling  pressures   predicted   by 

theories   reviewed   in    Section    I  are  plotted   versus   a  pipe   flexibility 

parameter 


EI 
.3 


in    Figures    7  and  8    for   loose  and   dense   soil   support, 
R  ' 
respectively.      For  comparison   purposes,    some  variables   have   been   held 

constant    in   all   cases.      A  smooth   steel    conduit    (v  =    0.3,    E  =    30   x    LO     psi, 

f      =    36,000  psi)   with  a   diameter  of    10    feet    is   assumed   to   be  covered   by 

10   feet   of   dry   sand  with   V     =0.4.      Loose   sand   supporting   the   culvert    is 

given  an   elastic   Young's   modulus   of   2000   psi,    while    the   modulus    for 

3 
dense   sand   is    taken    to   be  20,000  psi.      The  values   of   El/R     plotted   in 

Figures    7   and  8    reflect   conduit    stiffnesses    in   the    "flexible"    range, 

according   to    the    flexibility   parameter,    F    ,    given   by   Peck,    et   al    (1972), 

E     R"^    ,,  _   2, 


and  were  chosen  with  due  regard  to  commonly  manufactured  corrugated 
metal  sections. 

From  Figures  7  and  8  four  important  deductions  can  be  made: 

1.  All  theories  are  linear  on  a  log-log  plot. 

2.  The  mere  existence  of  surrounding  soil  increases  the 
critical  pressures  by  an  order  of  magnitude,  or  more, 
even  in  the  case  of  weak  supporting  soil. 
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3.  The  critical  buckling  pressure  is  strongly  influenced  by  the 
soil's  elastic  modulus.   If  this  modulus  is  increased  ten 
times,  the  critical  pressures  are  increased  by  roughly  the 
same  order  of  magnitude.   Thus,  cliaracterizing  the  quality  of 
backfill  by  percent  compaction,  or  by  rough  estimates  of  the 
soil  modulus,  is  patently  inadequate  from  the  standpoint  of 
predicting  the  critical  buckling  load. 

4.  Theory  indicates  that  for  weak  soil  support  buckling  is  the 
dominant  failure  mode  (compared  with  wall  crushing);  even 
for  strong  soil  support,  buckling  may  be  critical  in  the  more 
flexible  culverts.   Due  to  geometric  and  material  imperfections, 
and  to  residual  stresses  from  cold  forming,  observed  buckling 
pressures  may  be  even  lower  than  the  critical  pressures 
predicted  theoretically.   Accordingly,  buried  culverts  may  be 
more  susceptible  to  failure  by  buckling  than  is  now  generally 
acknowledged. 
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III.   BUCKLING  EXPERIMENTS  ON  SOIL  SUPPORTED  CONDUITS 

A.   Controlled  Experiments 
Because  of  the  complexity  of  the  interaction  between  flexible 
conduits  and  the  surrounding  soil,  and  the  vagaries  of  field  construction 
practices,  a  number  of  controlled  experiments  in  prefabricated  soil  bins 
have  been  conducted  to  compare  theory  with  observation.   As  in  any 
simulation  of  field  behavior  it  is  imperative  that  the  test  conditions 
closely  model  those  in  the  field.   Accordingly,  in  presenting  the 
results  of  a  number  of  investigations,  care  has  been  taken  to  review 
critically  the  apparatus  and  procedures  used  in  the  experiments. 

Bulson  (1963) 

Bulson   conducted   laboratory   tests   to    study    the   behavior  of   soil- 
supported,    flexible,    tin-plated   steel   conduits    (f     =   26   ksi).      Tubes   with  a 

thickness   of   0.011    inches  and   diameters   ranging   from  5   to   10   inches 

F  T 
{0.026   <  -^  <   0.208}  were   buried   in   a    five    foot    long   steel   test   cell 

R 

with  backfill  sand  compacted  by  a  vibrating  plate  in  three  inch  lifts. 

The  12  inch  long  tubes  were  restrained  at  the  ends  by  relatively  rigid 

observation  tubes  of  the  same  diameter  extending  axially  from  the 

flexible  tube  to  the  cell  wall.   Because  the  ratio  of  the  flexible 

length  to  the  diameter  was  not  sufficiently  large,  the  restraint 

significantly  increased  buckling  loads  over  that  expected  for  the 

plane  strain  case.   Tube  stresses  were  not  measured,  and  only  rough 

measurements  of  the  buckled  deformations  were  made. 

A  uniform  hydrostatic  pressure  was  applied  to  the  surface  of  the 

soil  cover,  which  was  varied  in  thickness  to  test  the  effect  of  burial 
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depth  on  conduit  performance.   The  ratio  of  burial  depth  (d)  to  the 

JO       o 

tube  diameter  (D)  was  either  "fT  ~  o"  or  7-»  and  the  corresponding  pressures 

producing  collapse  in  tubes  with  identical  geometry  were  compared. 

Surface  pressures  were  recorded  at  the  initiation  of  roof  collapse 

after  buckling  had  weakened  its  ability  to  carry  the  applied 

load.   Figure  9  illustrates  the  sequence  of  events  generally  observed 

during  loading  to  the  collapse  pressure,  p    •   Replication  of  p 

ma  X  ma  x 

was  difficult  and  average  values  were  used  in  comparisons  with  theoretical 
buckling  pressures  predicted  by  eqn  (36),  which  was  derived  by  Timoshenko 
and  Gere  (1961)  for  elastic  buckling  of  unsupported  cylinders  with 
clamped  ends  subjected  to  uniform  radial  pressure. 


Et        .     Et 
P   " 


""        R(n^~l){l+ 


nL 


2-,      3    "^ 
}   12R  (1-v  ) 


(,2_i)  ^  2n2-l-v 


ttR 


(36) 


The  theoretical  buckling  modes  (which  depend  upon  the  ratios  —  and  — )  were 

subsituted  into  eqn  (36)  to  calculate  p  .   It  was  found  that,  for  a  given 

J                             P  P 

d      ^.            ^  ■        max  J-    .    T                                  „,               .         max      ^          .                    ,       .    , 

— ,    the   ratio  was    fairly   constant.      The   ratio  also    increased  with 

Dp  p 

c  c 
d 
increasing  —  and   ranged   from   6.2    to    9.2.      Bulson  believed   these   results 

reflected   the   dependence  of   soil   support   on    the   depth  of    soil   cover. 

However,    because  of   the  nature  of   the  applied   load    (uniform  over   the 

entire   soil   surface),    the   soil   support   could  not   vary   much  over   the 

range  of   cover  heights    investigated.       Increased   soil  arching 

must   also   have   resulted   from   increased   soil   cover.      Thus,    increase 

of   the   collapse  overpressure  with   increasing   soil   cover  was  as 

indicative   of   a   change   in    the   load   distribution   on    the   tubes     as    it 

was   of   an   increase   in    the  modulus   of   soil   support. 
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(2)  Roof   Deflects 


(3)  Local  Buckles 


(4)    Buckle  Fails 


(5)  Roof  Collapses 


(6)    Rnal  Collapse 


FIGURE  9  STAGES  OF  TUBE  COLLAPSE  (after    BULS0N,I963) 
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Comparison  of   p  to   p      is    inappropriate.      Collapse   can   occur 

'  '^max  c 

immediately   after   buckling   is    initiated  or  at    loads    larger   than    tlie 

initial    buckling   load    (depending  on   soil   support   and   tube   stiffness). 

Comparison   of   the   initial   buckling  load    (rather   than   p        )   with   p 
^  max  c 

would  have   been   a   better   indicator   of    the  relation   between   the 
calculated  critical   pressure  and   the   observed   buckling   load.      Nevertheless, 
Bulson's   tests   clearly   demonstrated   the   importance   of   soil    support   on 
the   collapse  pressure. 

Luscher    (1966) 

Luscher  attempted  to  check  the  validity  of  previously  developed 
theoretical  buckling  solutions  by  modeling  experimentally  the  boundary 
conditions  and  types  of  loads  used  in  the  theoretical  analyses.   Two- 
ply,  bonded,  aluminum  foil  tubes  (f   estimated  to    equal  16  to  24  ksi) 
1.63  inches  in  diameter  and  6  inches  long,  with  a  composite  wall 
thickness  of  0.003  inches,  were  tested.   Each  smooth  tube  was  surrounded 
by  a  cylindrical  ring  of  Ottawa  sand,  either  3/8  in  or  2/3  in  thick. 
The  saturated  sand,  tested  in  both  loose  and  dense  states,  was  enclosed 
by  a  flexible,  impermeable  membrane  and  loaded  by  a  uniform  radial 
external  pressure  applied  to  the  membrane.   The  external  pressure 
producing  buckling  in  the  tube  was  compared  to  the  critical  pressure 
predicted  by  eqn  (8)  (Forrestal  and  Herrmann,  1965),  and  it  was 
foimd  that  the  theoretically  predicted  critical  pressure  was  much 
higher  than  that  measured  in  the  lab. 

Luscher  concluded  that  the  soil  support  modeled  by  Forrestal  and 
Herrmann  was  the  source  of  the  discrepancy,  and  consequently  he 
developed  an  empirical  expression  for  the  soil  support  modulus,  given 
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by  eqns  (11)  and  (12).   Radial  interface  pressures  acting  on  the  tube 
were  calculated  from 


P  =  P'  +  ^s  R  (3^) 


where,  p  =  local  radial  interface  pressure 

p'  =  uniform  external  pressure  applied  to  the  outside  of 
surrounding  soil  cylinder 

w  =  radial  tube  displacement  (positive  outwards) 
No  instrumentation  was  provided  to  measure  the  interface  pressures, 
thus  there  was  no  measured  verification  of  eqn  (37),   However,  a  check 
on  the  validity  of  eqn  (11)  was  made  by  back-calculating  K   from  each 
test  using  eqn  (10)  with  the  known  tube  geometry  and  failure  pressure. 
Luscher's  semi-empirical  theory  was  predicated  on  buckling  occurring 
in  high  modes,  but  actual  failure  configurations  seemed  to  take  the 
form  of  local  buckles.   Consequently,  the  mode  number  was  determined 
by  dividing  the  length  of  the  local  buckle  (I)    into  the  total  tube 
circumference.   Values  of  M  were  obtained  from  K  by  eqn  (11)  and 
compared  to  M  values  determined  from  oedometer  tests  on  Ottawa  sand 
at  the  densities  and  confining  pressures  used  in  the  buckling  tests. 
Good  correlation  was  found,  partly  because  the  formulation  of  eqn  (11) 
was  derived  from  the  test  observations.   No  appreciable  difference  in 
the  back-figured  support  modulus  was  found  between  loose  and  dense 
sand,  which  indicated  that  the  test  conditions  were  possibly  forcing 
artificial  behavior.   Luscher  believed  that  the  use  of  his  empirical 
soil  moduli  would  improve  the  predictive  capabilities  of  available 
theories  by  removing  a  large  amount  of  the  unconservative  error. 


55. 


However,  it  is  possible  that  tube  imperfections,  residual  stresses, 
and  non-composite  behavior  of  the  two-ply  aluminum  foil,  rather  than 
an  improper  formulation  of  soil  support,  were  causing  experimental 
failures  in  the  form  of  local  buckles  at  pressures  much  lower  than 
those  predicted  theoretically  for  buckling  in  high  modes.   The 
empirically  obtained  eqn  (11)  was  verified  only  for  the  tested  small- 
scale  geometries.   Luscher  admitted  that  scale  effects  could  be 
significant.   However,  he  believed  the  results  indicated  that  elastic 
buckling  was  an  important  potential  failure  mode  of  flexible  conduits 
supported  by  loose  backfill. 

Representative  results  obtained  by  Luscher  are  presented  in  Table  3. 
Replication  was  difficult,  as  illustrated  by  varying  buckling  pressures 
of  tubes  with  identical  stiffness  supported  by  sand  of  equal  density. 
Large  variation  in  the  observed  buckling  mode  was  probably  due  to 
differences  in  energy  released  during  snap-through  buckling,  which 
directly  affects  the  length  of  a  buckle  and  is  partly  dependent  on  the 
level  of  tube  imperfection.   Thus,  it  is  likely  that  local  imperfections 
in  the  tubes  tested  governed  their  behavior,  and  the  use  of  empirical 
relationships  obtained  from  these  tests  may  not  be  valid  for  field 
conditions. 

Allgood  and  Ciani  (1968) 

Allgood  and  Ciani  conducted  experimental  buckling  tests  on  smooth 
steel  (f  =  89  ksi)  tubes  buried  under  6  inches  of  uniform,  fine  sand. 

y 

The   tubes,    5    inches    in   diameter  and   20   inches   long  with  wall    thicknesses 
of   0.006,    0.012   and   0.018    inches,    were   tested   in   a   cylindrical    tank  with 
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a  diameter  of  5  Teet.   The  axis  of  the  tank  was  vertical.   The  walls 
were  constructed  to  minimize  friction,  and  the  base  of  the  tank  was 
sufficiently  removed  from  the  buried  tube  to  model  soil  of  infinite 
depth.   Soil  placement  with  density  control  was  accomplished  by  a 
raining  technique.   Load  was  applied  to  the  conduit  by  pressurizing 
a  pneumatic  bag  which  reacted  against  the  soil  surface. 

The  tubes  were  instrumented  with  electric  resistance  strain  gages 
placed  around  the  inner  and  outer  tube  circumference  to  measure  thrusts 
and  moments  that  developed  in  the  tube  walls.   A  rotating  deflection 
gage  was  installed  to  measure  the  inside  profile  of  the  tube,  and 
deflection  gages  were  placed  in  the  soil  below  the  tube  to  measure 
average  soil  strains. 

Tests  were  performed  with  varying  soil  densities  and  tube 
stiffnesses,  and  the  surface  soil  pressure  was  recorded  at  the  onset  of 
buckling.   In  all  tests,  strain  in  the  tubes  remained  elastic  prior  to 
buckling  and  failures  took  the  form  of  local  buckles  causing  almost 
immediate  collapse  in  thinner  tubes  and  collapse  after  increased  load 
in  thicker  tubes.   Interface  soil  pressures  were  calculated  in  two 
ways.   One  method  used  strain  measurements  to  determine  the  radius  of 
curvature  and  thrust  of  the  deformed  tube;  the  interface  pressure  was 
calculated  as  the  tube  thrust  divided  by  the  radius  of  curvature, 
assuming  small  bending  resistance  of  the  tube.   The  second  way  was  to 
calculate  the  interface  pressures  from  thrusts  and  moments  measured  in 
the  tubes.   The  pressures  calculated  from  these  methods  showed 
reasonable  agreement  before  buckling.   In  cases  where  the  tubes  buckled 
after  negligible  deflection,  pressures  were  larger  at  the  crown  than 
at  the  springline.   No  interface  pressures  were  measured  to  check  the 
accuracy  of  the  calculations.   Table  4  is  a  summary  of  the  experimental 
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results.   Replication  was  much  better  than  in  Luscher's  tests  but,  again, 
the  buckling  pressures  were  not  entirely  consistent  with  the  sand  density. 

Allgood  and  Ciani  related  the  interface  buckling  pressure  to  the 
surface  soil  pressure  by  an  arching  factor, 

p^  =  P3  (1  -  C)  (38) 

p  =  applied  surface  soil  pressure 

C  =  arching  factor  =  {l  -  F  /p  R},  where  F   is  the  measured  thrust  at 

c   s  c 

the  springline. 
Arching  was  found  to  vary  among  tests  with  identical  soil  densities 
and  tube  stiffnesses  and  this  was  thought  to  be  due  to  improper  seating 
of  the  tube  during  placement.   The  calculated  interface  pressure  was 
compared  to  the  critical  uniform  external  buckling  pressure  predicted  by 
Chelapati  and  Allgood  {eqn  (18i)},  and  fair  agreement  was  found,  although 
eqn  (18i)  was  developed  for  high  mode  buckling  and  only  local  buckling 
was  observed. 

In  the  interpretation  of  buckling  tests  the  writers  foimd  it  useful 
to  compare  the  buckling  stress,  f  ,  to  the  yield  stress  of  the  wall 
material,  f  .   An  efficiency  ratio,  R  ,  was  defined 

f 
Rf  =  ^  (39) 

y 

as  a  measure  of  how  efficiently  the  potential  load  capacity  of  the  conduit 
was  being  utilized.   In  cases  where  f  was  not  measured  in  the  vicinity  of 
the  buckle   it  was  calculated  using  the  best  technique  available, 
consistent  with  the  nature  of  the  input  data.   The  values  of  R  in 
Table  4  clearly  show  that  elastic  buckling  occurred  at  pressures  well 
below  those  required  to  induce  failure  by  wall  crushing. 
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From  the  measured  behavior  of  the  loaded  tubes  Allgood  and  Cianl 
concluded  that: 

(1)  deflections  of  buried  flexible  cylinders  under  a  given 
load  are  governed  mainly  by  the  soil  modulus, 

(2)  interface  soil  pressures  can  be  determined  from 
measurements  of  strains  in  the  cylinder,  and 

(3)  buckling  loads  may  be  predicted  using  eqn  (18i). 

The  modulus  of  elastic  support  used  in  eqn  (18i)  was  determined 
using  eqn  (11)  developed  from  Luscher's  buckling  experiments.   The  use 
of  these  equations  allowed  correlation  between  theory  and  experiment, 
but  the  critical  interface  pressure  assumed  in  eqn  (18i)  was  uniform 
while  the  interface  pressures  calculated  from  measured  strains  were 
not.   In  addition,  failure  by  elastic  buckling  in  high  modes,  modeled 
by  eqn  (ISi),  was  not  observed  in  the  experiments.   These  factors, 
combined  with  the  small  scale  of  the  experiments,  make  it  unlikely 
that  the  test  results  can  be  scaled  to  larger  conduits. 

Dorris  (1965)  and  Albritton  (1968) 

Studies  at  the  Waterways  Experiment  Station  by  Dorris  and  Albritton 
investigated  the  effects  of  burial  depth  on  the  buckling  behavior  of 
smooth,  flexible  tubes.   The  tubes  were  instrumented  so  that  quantitiative 
data  could  be  obtained  to  supplement  Bulson's  (1963)  observations.   Dry 
sand  at  a  relative  density  equal  to  80%  was  used  to  cover  the  tubes  with 
depths  from  0-2  diameters.   Aluminum  cylinders  with  diameters  of  3.5 
inches,  lengths  10.5  and  12  inches,  and  wall  thicknesses  of  0.01,  0.022 
and  0.065  inches  were  tested  in  a  cell  with  special  provisions  made  to 
reduce  soil  friction  along  the  cell  wall.   The  cell  length  was  4  feet 
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and  the  base  of  the  cell  was  2  feet  below  the  cylinder.   The  cylinder  ends 
were  enclosed  by  a  flexible  gasket  and  an  aliiminum  cap  so  that  plane 
strain  conditions  were  approximated. 

Instrumentation  included  metal  foil  strain  gages  installed  inside 
and  outside  of  the  tube  at  the  crown,  invert,  and  springline.   This 
allowed  determination  of  the  thrusts  and  moments  in  the  tube  at  the 
points  of  instrumentation,  and  these  were  compared  at  different 
tube  burial  depths.   Uniform  pressures  were  applied  to  the  soil  surface. 
The  recorded  critical  pressure  was  the  surface  overpressure  at  which 
a  buckle  first  formed  in  the  tube. 

Failure  modes  were  observed  to  vary  with  burial  depth,  partially 
agreeing  with  Bulson's  test  results.   For  burial  of  0  to  y  diam  (D), 

failure  took  the  form  of  a  catastrophic  snap-through  buckle  at  the 

3 
crown.   Burial  at  t-  D  changed  the  location  of  the  local  buckle  to  the 

springline,  and  tubes  buried  at  greater  depths  buckled  locally  at  the 

invert.   Measured  strains  before  buckling  were  elastic  in  all  cases, 

and  the  strains  measured  at  the  initiation  of  buckling  varied  with  the 

depth  of  soil  cover.   The  largest  strains  were  measured  at  the  crown  when 

shallow  cover  (less  than  —   D)  was  provided.   With  increasing  depth  of 

cover,  the  largest  measured  strains  were  transferred  from  the  crown 

to  the  springline,  and  then  to  the  invert.   Thus,  strains  measured  in 

these  tests  verified  that  variation  of  soil  cover  significantly  affects 

the  load  distribution  on  conduits  with  shallow  (less  than  one  diameter)  burial. 

Results  from  a  few  tests  conducted  by  Dorris  using  compacted 
(approximately  95%  Std.  AASHO)   highly  plastic  clay  as  backfill  showed 
that  eqn  (2)  from  Timoshenko  and  Gere  (1961)  was  not  a  bad 
approximation  for  the  buckling  stress,  indicating  that  either  the 
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backfill  provided  little  buckling  restraint  or  that  the  small  tubes 
possessed  imperfections  that  grossly  affected  the  test  results. 

Large  scatter  was  reported  in  the  experimental  results,  due 
partly  to  instrumentation  errors.   Dorris  also  found  that  the  sand 
placement  technique  (raining  vs.  rodding)  influenced  the  test  results. 
Sand  placed  at  a  given  density  by  raining  provided  less  support  than 
placement  by  rodding.   This  could  have  been  due  to  the  inability  of 
the  raining  technique  to  provide  the  required  density  immediately 
around  the  tube,  and  to  the  pre-straining  effects  produced  by  rodding, 

Watkins  and  Moser  (1969) 

Watkins  and  Moser  reported  the  results  of  controlled  tests  on  3 
to  5  foot  diameter,  and  20  foot  long  steel  culverts  (f  =  40  ksi) 
with  2  2/3  X  1/2  inch  and  3x1  inch  corrugations  and  16  and  18  gage 
wall  thicknesses. 

A  test  bin  24  feet  long,  15  feet  wide  and  18  feet  high  was  used 
to  contain  the  culvert  and  backfill,  and  at  least  five  feet  of 
clearance  on  all  sides  of  the  conduit  was  provided.   The  bin  walls 
were  elliptical  to  reduce  side  friction  effects.   Dry  silty  sand 
compacted  to  densities  ranging  from  61  -  103%  Std.  AASHO 

(7  5  to  129  pcf)  was  used  as  backfill. 

Loads  were  applied  by  hydraulic  jacks  bearing  against  a  steel 
plate  on  the  surface  of  the  soil  5  feet  above  the  culvert  crown. 
Culverts  were  instrumented  to  measure  deformations  in  cross-section. 
Several  pressure  gages  were  also  placed  around  the  circumference  to 
measure  normal  pressures  near  the  soil-culvert  interface. 
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The  purpose  of  the  investigation  was  to  establish  performance 
limits  for  buried  culverts  and  to  develop  a  means  for  estimating  when 
a  performance  limit  will  be  reached,   A  summary  of  the  published 
experimental  results  is  given  in  Table  5. 

From  their  test  results  and  observations  Watkins  and  Moser  concluded 
that: 

(1)  The  performance  limit  should  be  defined  as  the  deformation 
beyond  which  either  the  culvert  or  soil  can  no  longer 
perform  its  design  function.   Usually  this  occurs  when  the 
resistance  to  additional  deflection  approaches  zero. 

(2)  In  dense  soil  the  performance  limit  was  determined  by  the 
crushing  strength  of  the  culvert  wall. 

(3)  Buried  culverts  are  capable  of  sustaining  larger  loads 
after  local  buckling  is  initiated. 

Empirical  curves  were  developed  relating  the  thrust  stress  in 
the  culvert  wall,  calculated  by  eqn  (3)  at  the  interpreted  performance 
limit,  to  the  slendemess  ratio  of  the  pipe  wall  for  different  soil 
support  conditions  (in  terms  of  percent  Std.  AASHO  maximum  density). 
The  curves  were  designed  to  reflect  the  fact  that,  with  good  soil 
support,  as  the  slendemess  ratio  decreased  the  performance  limit 
passed  successively  from  wall  crushing,  through  inelastic  buckling,  to 
elastic  buckling.   In  some  cases  the  thrust  stresses  at  failure 
calculated  from  eqn  (3)  exceeded  the  yield  stress.   It  was  decided 
that  this  was  due  to  neglect  of  soil  arching,  and  load  factors  were 
developed  to  adjust  the  calculated  stresses  to  the  yield  stress. 
These  load  factors  were  related  solely  to  the  relative  compaction  of 
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the  soil;  the  influence  of  the  conduit  stiffness  was  ignored.   No 
measurements  were  made  to  verify  the  calculated  stresses  or  to 
differentiate  between  inelastic  and  elastic  buckling.   Performance 
limits  were  based  primarily  on  visual  interpretations  which,  as 
discussed  earlier,  can  be  misleading.   The  above  mentioned  empirical 
curves  developed  from  the  test  results  have  been  adapted  for  use  in 
design  practice  (Handbook  of  Steel  Drainage  and  Highway  Construction 
Products,  1971). 

Meyerhof  and  Baikie  (1963) 

Meyerhof  and  Baikie  (1963)  tested  segments  of  steel  cylinders 
bearing  against  dense  sand  to  verify  their  formula  for  the  critical 
elastic  buckling  stress  of  buried  flexible  culverts  (eqn  2  7).   Smooth 
and  corrugated  (3/4  x  1/4)  quadrantal  sections  with  radii  of  12  to 
24  inches   and  wall  thicknesses  of  0.017  and  0.042  inches   were  loaded 
with  a  compressive  thrust  on  one  of  the  longitudinal  edges,  while  the 
other  longitudinal  edge  was  held  fixed.   Electric  resistance  strain 
gages,  soil  pressure  gages,  and  deflectometers  were  mounted  along  the 
central  circumference  of  the  plate.   The  tests  were  conducted  in  a 
steel  bin  two  feet  long,  four  feet  wide,  and  three  feet  deep.   Since 
the  axial  length  of  the  quadrantal  plate  was  19  inches,  only  2  1/2 
inches  of  sand  separated  the  ends  of  the  sections  from  the  bin.   At 
least  one  diameter  of  sand  backfill  compacted  to  90%  relative  density 
and  surcharged  v/ith  approximately  2  psi  overpressure  separated  the 
plate  section  laterally  from  the  bin  wall.   No  special  effort  was  made 
to  reduce  friction  between  the  sand  and  the  bin,  hence  the  boundary 
effects  of  the  bin  were  not  eliminated.   Some  test  results  obtained  by 
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Meyerhof  and  Baikie  are  summarized  in  Table  6.   Replication  was  difficult, 
especially  in  tests  on  the  most  flexible  conduit  in  which  the  range  in 
buckling  stress  was  approximately  70%  of  the  mean. 

Stresses  calculated  from  measured  strains  in  the  loaded  plates 
showed  that  the  more  flexible  smooth  plates  buckled  at  approximately 
15%  of  the  yield  stress,  and  the  least  flexible  corrugated  plates  failed 
by  crushing  near  the  yield  stress.   Average  radial  strains  before 
buckling  were  less  than  1%.   Soil  pressure  measurements  indicated  that  the 
normal  interface  soil  pressures  were  not  uniform  prior  to  initiation  of 
buckling. 

Although  it  was  claimed  that  eqn  (2  7)  accurately  predicted  the 
experimental  results, 

(1)  the  modulus  of  soil  restraint  was  not  constant  along  the 
plate-soil  interface  at  the  initiation  of  buckling,  as 
assumed  in  the  theory;  and 

(2)  high  modes  of  buckling  were  not  observed  in  the  tests. 
Instead,  failure  always  occurred  in  the  form  of  a  local 
buckle  near  the  soil  surface  because  of  lack  of  soil 
restraint  inherent  in  the  experimental  apparatus. 

In  tests  on  corrugated  plates  with  constant  wall  stiffness  (EI) 
the  efficiency  ratio  increased  as  the  radius  of  curvature  decreased, 
which  is  the  inverse  of  the  behavior  observed  by  Watkins  and  Moser 
(see,  for  example,  the  second  and  fourth  tests  listed  in  Table  5,  in 
which  R^  increased  with  increasing  radius).   This  inconsistency  can  be 
related  to  the  hinged  end^condltion  which  did  not  allow  full 
development  of  bending  moments  in  the  zone  of  low  soil  restraint  where 
buckling  always  occurred.   This  limited  the  failure  mode  either  to  wall 
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crushing  or  elastic  buckling;  thus,  as  the  plate  radius  decreased, 
elastic  buckling  resistance  increased  and  a  larger  efficiency  ratio 
resulted.   In  Watkins  and  Moser's  tests  bending  moments  were  allowed 
to  develop.   Because  decreasing  the  culvert  radius  tends  to  increase 
the  bending  stress  relative  to  the  thrust  stress,  it  is  likely  that 
buckling  was  initiated  at  loads  below  those  required  to  cause  buckling 
in  the  absence  of  bending  stresses.   Accordingly,  the  load  carrying 
efficiency  (reflected  by  the  efficiency  ratio)  was  reduced.   This 
important  aspect  of  conduit  behavior  was  not  reflected  by  Meyerhof 
and  Baikie's  tests. 

\ 
Howard  (1971) 

Howard  conducted  tests  on  uncorrugated,  7  foot  long  steel  conduits 

supported  by  sandy  clay  compacted  to  90%  and  100%  Std.  AASHO  density. 

The  diameters  tested  ranged  from  18  to  30  inches,  and  the  pipe 

EI 
flexibility  parameter  — -■  varied  from  0.30  to  10.2  psi.   The  steel 

R 

soil  container  within  which  the  conduits  were  buried  was  6  feet  wide, 

7  feet  long  and  7  feet  deep.   The  ends  of  the  conduits  extended  through 

the  container  walls,  and  provisions  were  made  to  prevent  appreciable 

restraint  to  deflection.   The  container  walls  were  covered  with  a 

polyethylene  liner  lubricated  with  grease  to  reduce  interface  friction. 

The  longitudinal  axis  of  each  conduit  was  placed  four  feet  below  the 

soil  surface.   A  thick  wooden  load-plate  was  placed  on  the  soil  surface 

to  distribute  the  surcharge  load  applied  by  a  large  universal  testing 

machine.   Instrumentation  on  each  conduit  included  SR-4  strain  gages 

measuring  strains  on  the  inside  surface,  a  rotating  dial  deflectometer 

measuring  the  deflected  profile,  and  soil  pressure  cells   installed 

flush  with  the  outside  surface  at  the  crown,  invert  and  springlines. 
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measuring  normal  interface  pressures.   Soil  displacement  indicators 
were  also  positioned  in  the  backfill. 

Test  results  are  given  in  Table  7.   Observed  behavior  during  loading 
varied,  depending  on  the  soil  density  and  conduit  stiffness.   In  the 
weaker  soil  (90%  Std.  AASHO  density)  deflections  were  generally  large 
under  relatively  small  surcharge  loads.   The  most  flexible  conduit 
experienced  two  local  elastic  buckles  after  8%  deflection,  which  marked 
the  performance  limit.   The  less  flexible  conduits  formed  small  inelastic 
crimps  after  measured  strain  on  the  interior  of  the  conduit  wall  exceeded 
the  yield  strain.   A  performance  limit  was  reached  after  larger  surcharge 
loads  were  applied  and  was  characterized  by  the  formation  of  a  large 
inelastic  snap-through  buckle,  usually  at  the  crown.   Because  of  the 
instrumentation  scheme,  it  was  impossible  to  determine  the  magnitude  of 
developed  moments,  hence  their  Influence  on  the  failure  mode 
could  not  be  evaluated.   Generally,  as  the  conduit  stiffness  increased 
the  number  of  observed  inelastic  crimps  decreased  and  the  load  at  which 
they  formed  increased.   In  all  cases  when  weak  backfill  was  provided, 
inelastic  crimping  occurred  outwards  into  the  backfill. 

With  the  denser  backfill  (100%  Std.  AASHO  density)  deflections  were 
reduced.   However,  local  elastic  buckling  occurred  in  the  more  flexible 

conduits  at  approximately  the  same  loads  which  produced  snap-through 

EI 
buckling  in  weaker  backfill.   Only  in  the  less  flexible  conduits  (— r-  >  4) 

R 

did  the  improved  backfill  prevent  wall  buckling  before  the  capacity  of  the 

apparatus  was  exceeded.   Howard  concluded  that  although  denser  backfill 

reduced  deflection,  it  did  not  significantly  increase  the  critical 

buckling  load.   However,  the  writers  point  out  that  the  load  capacity  of 

the  conduit  supported  by  the  denser  backfill  is  well  above  that  of  the 

same  conduit  supported  by  weaker  soil. 
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Measurement  of  the  normal  interface  soil  pressures  at  the  initiation 
of  buckling  showed  that  in  most  cases  the  pressure  distribution  was  not 
uniform.   When  appreciable  deformation  of  the  conduit  occurred  before 
buckling,  the  crown  pressure  was  up  to  40%  less  than  the  springline 
pressure.   If  buckling  occurred  after  relatively  small  deformation  the 
pressure  distribution  was  more  uniform.   Some  of  the  pressure  measurements 
were  unrealistic,  illustrating  the  difficulty  of  obtaining  good  earth 
pressure  measurements,  especially  near  the  soil-conduit  interface. 

Soil  displacement  indicators  showed  that  most  displacement  occurred 
within  one  radius  of  the  conduit  wall.   Although  almost  three  radii  of 
soil  separated  the  container  walls  from  the  conduit,  it  was  noted  that 
the  container  did  influence  soil  displacements.   Howard  believed  that 
the  effects  were  not  large  enough  to  control  behavior,  but  he  suggested 
that  field  tests  were  needed  to  correlate  and  refine  empirical  relation- 
ships developed  from  laboratory  tests.   Nevertheless,  Howard's  tests 
demonstrated  that  even  well-compacted  clay  backfills  (100%  Std.  AASHO) 
provide  inferior  support  to  that  of  granular  soils,  especially  from  the 
standpoint  of  controlling  deflections. 

B.   Field  Experiments 
Verification  of  theoretical  models  predicting  conduit  behavior, 
especially  failure  modes,  must  involve  comparisons  with  field  performance. 
The  number  of  field  tests  that  can  be  used  for  this  purpose  is  limited 
because  instrumentation  is  seldom  provided  to  measure  all  of  the  soil 
properties  and  conduit  loads  needed  as  input  for  most  predictive  models. 
Moreover,  field  installations  are  costly  and  hence  are  seldom  designed 
to  reach  failure.   Thus,  the  failure  mechanisms  are  unclear  and  the 
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factors  leading  to  failure  are  uncertain.   The  investigations  reviewed 
below  have  been  instrumented  and  documented  sufficiently  to  provide 
insight  into  the  field  behavior  of  buried  conduits. 

Prairie  Farm  Rehabilitation  Administration  (1957) 

To  observe  the  behavior  of  flexible  culverts  in  clay  backfill,  a 

34  feet  high  test  fill  covering  a  6  foot   diameter  steel  culvert  (f  =  40  ksi) 

FT 
was  constructed.   The  20  gage,  (2  2/3  x  1/2)  corrugated  culvert  (~t  =  0.7} 

R 

was  instrumented  with  strain  gages  to  measure  both  compressive  and  bending 

strain  and  a  field  protractor  was  used  to  measure  the  deflected  shape. 

A  clay  backfill  with  medium  plasticity  was  compacted  to  98%  Std.  AASHO 

density.   The  com.pacted  soil  unit  weight  around  the  culvert  showed  +  11% 

variation  about  the  mean,  indicating  nonuniform  soil  support.   After  the 

final  fill  height  was  achieved,  a  small  outward  buckle  was  observed  in 

the  lower  haunch.   The  average  measured  thrust  stress  was  20,000  psi 

(R  =  0.5),  and  bending  stresses  were  sufficiently  large  to  initiate 

outer  fiber  yielding  at  the  location  of  the  buckle.   A  vertical 

deflection  of  4  inches  (5.5%  of  the  diameter)  was  measured  just  before 

the  formation  of  the  buckle.   It  was  speculated  that  poorly  compacted 

soil  at  the  haunch  provided  a  zone  of  weak  support  which  allowed  local 

buckling.   Despite  the  formation  of  a  local  buckle,  the  culvert  had  not 

yet  reached  its  load  capacity. 

A  second,  similarly  instrumented  culvert  with  a  9  foot  diameter, 

FT 
1  gage  wall  thickness  and  6x2  inch  corrugations  {-^  =    30.5}  was 

R 

backfilled  with  a  gravelly  clay  compacted  to  94%  Std.  AASHO  density  and 

loaded  with  a  99  foot  high  fill.   The  overburden  load  caused  a  vertical 

crown  deflection  measuring  3.5%  of  the  culvert  diameter  and,  combined 
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with  an  average  compressive  thrust  stress  of  12,500  psi,  produced 
moments  large  enough  to  initiate  yielding  at  several  locations, 
although  no  wall  disturbance  was  observed. 

Demmin  (1966) 

Demmin    reported   the   behavior  of  a   long-span   pipe   arch   subjected  to 
live   loads  applied   through   5.2    feet   of   soil   cover    (one-quarter  of   the 
span).      The  pipe  arch,    with  a   span   of   20   feet    7   in   and  a   rise  of   13   feet 

2    in.,    was   constructed  of   7   gage,    6x2    inch  corrugated  curved  plates 

I 
3 


FT 
{—  =  1.5}.   The  pipe  was  placed  within  an  excavation  having  very  firm 

R 
side  slopes,  and  backfill  placed  around  the  pipe  consisted  of  a  sandy 

gravel  compacted  to  107%  Std.  AASHO  density.   Instrumentation  included 

strain  gages  positioned  on  adjacent  corrugation  crests,  troughs,  and 

neutral  axes  at  six  locations  along  the  circumference,  which  enabled 

determination  of  the  thrusts  and  moments  in  the  pipe  wall  at  these 

locations.   Deflections  were  measured  with  phototheodolites,  and  three 

Heierli  earth  pressure  cells  were  installed  in  a  horizontal  plane  four 

inches  above  the  pipe  crown  to  measure  the  percentage  of  the  live 

load  transferred  to  the  pipe.   The  live  load  was  applied  by  piling 

steel  slabs  on  an  88  square  foot  grillage  of  railroad  ties  resting 

on  the  soil  surface  above  the  crown. 

Strains  produced  by  a  690  ton  surface  load  approached  yield  in  two 

locations  on  the  upper  half  of  the  pipe.   However,  residual  strains 

produced  during  placement  and  backfilling  of  the  pipe  combined  with  the 

live  load  strains  to  give  a  net  strain  well  within  the  elastic  limit. 

Bending  stress  comprised  65%  of  the  load-induced  extreme  fiber  stress, 

and  vertical  deflection  was  1.3%  of  the  span  (3.1  in).   t\Fhen  the  live 
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load  was  increased  to  850  tons,  two  inward  bulges  with  amplitudes  of 
11.8  and  5.9  inches  formed  on  either  side  of  the  crown  near  the  location 
of  the  maximum  moments.   Because  the  strain  gages  had  been  removed, 
the  state  of  stress  in  the  pipe  was  unknown.   The  vertical  deflection 
was  approximately  1.7%  of  the  span  (4  in).    Vertical  earth  pressures 
measured  at  the  crown  reached  65  psi.   After  the  surface  load  was 
increased  to  954  tons,  the  grillage  was  observed  to  settle  and  cause 
spreading  of  the  load  over  a  larger  surface  area.   This  relieved  the 
pressure  on  the  pipe  arch  and  caused  the  steadily  increasing  vertical 
deflection  to  stabilize  at  2.7%  of  the  span  (6.6  in).    It  was  believed 
by  Demmin  that  in  the  absence  of  spreading,  954  tons  would  have  been 
close  to  the  collapse  load. 

The  test  demonstrated  that: 

(1)  bulging  can  occur  in  long-span  culverts  due  to  the 
influence  of  bending  moments  at  thrust  stresses  well 
within  the  elastic  limit, 

(2)  very  strong  soil  support  can  prevent  buckling  and 
collapse  despite  the  localized  occurrence  of 

of  large  plastic  deformations, 

(3)  collapse  of  long-span  culverts  due  to  excessive 
deflection  may  occur  at  crown  deflections  below 
5%  of  the  span,  and 

(4)  stresses  induced  during  installation  of  long-span 
culverts  with  excellent  soil  support  can  be  large 
and  may  even  approach  the  yield  stress.   If 
compaction  is  symmetric  and  yielding  is  not 
excessive,  the  bending  stresses  produced  during 
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installation  are  opposite  to  those  induced  by  surcharge 
loading  and  may  have  a  favorable  effect  on  performance. 


Wu  (1977) 


Wu  examined  the  performance  of  several  culverts  that  were  installed 
without  modem  compaction  control.   A  device  designed  to  measure  in  situ 
soil  moduli  was  used  to  determine  the  level  and  uniformity  of  soil 
support  obtained  by  the  backfilling  operations.   Wide  variation  of  the 
interpreted  elastic  moduli  was  found  around  each  culvert.   Strain  gages 
were  placed  on  the  inside  of  the  culverts  and  small  coupons  containing 
the  gages  were  carefully  removed  from  the  wall.   The  measured  changes  in 
strain  were  assumed  to  be  those  induced  by  the  overburden  load  plus 
residual  strains  produced  during  fabrication.   Limited  measurements  of 
residual  strains  in  recently  fabricated,  4  foot  diameter  corrugated 
culverts  were  made,  but  the  results  showed  wide  scatter.   Even  so,  the 
importance  of  residual  stresses  was  made  evident  by  the  fact  that  their 
average  value  approached  8000  psi;  hence,  if  their  influence  is  not 
considered,  comparisons  between  theoretical  and  measured  behavior  are 
dubious. 
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IV.   DISCUSSION  OF  EXPERIMENTAL  RESULTS 

A.   Comparison  of  Results  from  Controlled  Buckling  Tests 
Results  of  controlled  experiments,  in  which  the  buckling  pressure 
and  mode  have  been  identified,  are  plotted  in  Fig.  10  in  terms  of  a 
scaled  buckling  pressure  vs.  a  soil -conduit  flexibility  factor. 
There  is  considerable  scatter,  but  the  trends  are  clear: 

(1)  Kloppel  and  Clock's  buckling  theory  follows  the  pattern 
of  test  results  tolerably  well;  scatter  is  greatest 
when  the  mode  of  failure  is  elastic  buckling. 

(2)  Imperfeccions  and  residual  stresses  in  the  tested 
conduits  lower  the  buckling  pressure,  particularly 
in  the  case  of  small  scale  tests  (e.g.  Luscher's 
experiments). 

(3)  The  use  of  high  yield  point  steel  (Allgood  and  Ciani's 
experiments)  increases  the  load  at  which  wall  crushing 
and  plastic  yielding  will  occur  and  may  make  buckling 
more  critical  even  with  excellent 

soil  support.   At  the  same  time,  the  effect   of 
imperfections  and  residual  stresses  on  the  buckling 
pressure  is  reduced.   Accordingly,  the  capacity  of 
the  conduit  to  carry  loads  is  increased  significantly; 
the  potential  for  savings  in  cost  and  material  resources 
needs  to  be  examined  by  full  scale  experimentation. 

B.   Important  Deductions  Regarding  Conduit  Behavior 
In  spite  of  the  studies  that  have  been  made  to-date,  there  remains  a 
large  gap  between  the  results  of  experiments  and  their  theoretical 


77. 


E°+  \ 

m 

•a 

1-1  + 

\ 

; 

0) 

u 

3 

m 

x: 
en 

3 

■ 

CO 

t-i 

+ 

\ 

s« 

•J 
CD 
O 

O 

i- 

Cl 
O 

u 

3 

m  1       1 

\ 

3 

c 

a 

+ 

s+ 

\ 

«• 

o 

'♦- 
in 
O 

Ul 

tn 
O 

LlI 
1 

'in 
O 

a> 

c 

1 

B+ 

\ 

-1 

tn 

a 

<o+ 

ew\ 

O 

(/^« 

N 

0} 

Ji 

^ 

o 

/ 

1 

y' 

o 

/^ 

/ 

/ 

o/ 

u 

X 

® 

q5 

(/)  \ 

Q. 

a 

CO  Q 

□\ 

:o 

. 

^ 

Ad 

-1 
► 

"► 

-1 

► 

\-j>V 

w<;J 

*"-* 

\ 

tn 

CO 

to 
(J) 

0^ 

to 

0} 

WQ      "'Q 

\ 

ir> 

3 

^ 

3 

^^ 

a» 

\ 

^^y 

0) 

q: 

U) 

a: 

-^ 

v)e 

\ 

V 

"1       ■ 

0} 

o 

O 

5 

CO 

(OlD 

\ 

j^ 

OD 

^ 

CD 

M— 

\ 

T3 
O 
O 

a 

c 

O 

x: 

1- 
o 

\ 

\ 

u 

o 

3 

C7> 
< 

1 

5 

1 

X 

\     ^> 

"5 

o> 

\ 

E 

< 

a 

o 

• 

► 

+ 

en  □ 

\ 

to 

_J 

< 

1 

\-J  _i  -i-J 

-  o 


LU 


CL 
X 

UJ 


o 

ct    CQ 

UJ     Q 

LJ 
-      _J 

:    o 


-  o     - 


u 


-  o 


o 


o 
o 


'3 


o 
o 
C5 


O 
O 


O 


C/) 

z> 
(/^ 

LU 

q: 

Ll 

o 

z 
o 
02 

2 

o 
o 


(/^ 


Q 

•z 
o 
o 

Q 
LU 

CC 
CO 


UJ 


78. 


interpretation,  which  limits  their  use  in  developing  generalized, 
rational  design  procedures.   The  shortcomings  of  available  test  results 
are  examined  below  with  a  view  towards  identifying  future  studies  most 
needed  for  the  development  of  improved  design  approaches. 

Modulus  of  Soil  Restraint 


Experimental  results  clearly  demonstrate  that  the  mode  of  failure, 
and  hence  the  performance  limit,  is  strongly  influenced  by  the  nature 
(linear,  nonlinear,  time -dependent)  and  magnitude  of  the  support  provided 
by  the  soil  backfill.   All  high  mode  buckling  theories  known  to  the 
writers  assume  the  soil  to  be  linear,  elastic  and  time- independent  in 
its  behavior.   The  analysis  of  snap-through  buckling  by  Klb'ppel  and 
Clock  considers  stress-dependent  soil  support,  but  the  modulus  is 
assumed  to  increase  linearly  with  the  stress  level.   Because  the  true 
nature  of  soil  support  has  not  yet  been  modeled  appropriately,  no 
buckling  theory  untempered  by  empirical  correlations  can  be  a  guide  to 
practical  design.   On  the  other  hand,  empirical  correlations  dealing 
with  the  soil  support  modulus  are  based  on  the  results  of  experiments  on 
small  tubes  (e.g.  Luscher,  1966;  Allgood  and  Ciani,  1968).   As  the 
soil  modulus  is  very  sensitive  to  the  type  of  soil  (all  the  above 
experiments  were  conducted  on  sands)  and  on  the  level  of  stress  and 
strain  in  the  soil,  the  results  of  tests  on  small  tubes  cannot  be 
scaled  directly  to  field  conditions.   Only  in  the  studies  by  Wu  (1977) 
was  there  any  attempt  to  measure  the  actual  soil  modulus  in  backfill 
supporting  large  culverts,  but  the  results  are  typical  only  of  culverts 
constructed  before  controlled  compaction  was  commonly  specified.   In 
future  experiments,  the  best  efforts  possible  should  be  directed 
towards  measuring  stress  and  strain  in  soil  backfill  surrounding  a 
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conduit.   Tests  should  also  be  conducted  using  clay  backfills  to 
investigate  their  support  characteristics  as  compared  to  granular 
backfills,  since  it  is  not  always  feasible  to  specify  the  use  of 
granular  soils  in  practice.   Experimental  measurements  of  normal 
interface  pressures  are  limited  in  number,  and  many  of  the  results 
obtained  are  suspect.   The  buckling  load  is  dependent  on 
the  distribution  of  boundary  pressures,  and  it  is  important  to  determine 
experimentally  what  distributions  may  be  expected  and  how  they  may 
influence  the  buckling  load.   Shear  stresses  at  the  soil-conduit 
interface  have  been  considered  only  approximately  in  theory  and 
largely  ignored  experimentally.   Efforts  should  be  made  to  clarify 
their  role  in  conduit  performance. 

Imperfections  and  Residual  Stresses 

It  is  not  possible  to  fabricate  conduits  without  some  material  and 
geometric  imperfections  and  without  leaving  residual  stresses  due  to 
forming  (unless  they  are  relieved  by  heat  treating).   Imperfections  and 
residual  stresses  influence  the  mode  of  failure  in  ways  not  directly 
scalable  from  tests  on  small  tubes.   They  complicate  efforts  to  obtain 
reproducible  experimental  results,  and  the  dangers  inherent  in  drawing 
conclusions  from  the  results  of  a  few  experiments  (let  alone  from  a 
single  field  experience)  cannot  be  emphasized  too  strongly. 

Strength  tests  performed  on  corrugated  (2  x  1/2  and  2  2/3  x  1/2) 
steel  culverts  by  Lane  (1965)  showed  the  influence  of  imperfections  on 
culvert  behavior.   Over  150  tests  on  culverts  ranging  from  12  to  48 
inches  in  diameter  and  from  0.052  to  0.168  inches  in  wall  thickness 
were  conducted  to  investigate  the  ability  of  corrugated  pipe  to  resist 
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yielding,  buckling,  and  seam  slippage.   The  test  apparatus,  illustrated 
in  Fig.  11,  encouraged  wall  disturbance  at  gaps  in  circumferential 
support  located  at  the  springlines.   The  buckling  load  was  observed  to 
increase  as  the  span  decreased,  until  a  limiting  size  was  reached  such 
that  the  bucklijig  load  was  no  longer  dependent  on  the  gap,  and  local 
buckling  occurred  randomly  around  the  pipe  circumference  instead  of  at 
the  springline.   Lane  interpreted  the  cause  of  this  behavior,  as  well 
as  the  scatter  in  experimental  results,  to  be  geometric  Imperfections  in 
the  pipe  x-sections.   The  writers  believe  that  reducing  the  gap  span  to 
the  limiting  size  provided  sufficient  circumferential  support  to  encourage 
high  mode  buckling, but  imperfections  in  the  form  of  geometric  irregularities 
or  residual  stresses  caused  local  buckling  to  occur  at  random  locations 
rather  than  in  a  large  number  of  circumferential  waves.   Because  the 
applied  buckling  load  could  not  be  compared  to  a  theoretical  critical 
load,  due  to  uncertain  load  distribution  in  the  tests,  the  reduction 
in  buckling  resistance  due  to  imperfections  can  not  be  determined.   In 
the  future,  tests  should  be  conducted  on  full-size  culverts  with  special 
efforts  directed  towards  evaluating  the  effects  of  imperfections  and 
residual  stresses  on  the  buckling  load. 

Bending  Stresses  and  Local  Yielding 

The  economy  inherent  in  the  use  of  flexible  metal  culverts  is 
predicated  on  sufficient  deflection  developing  so  that  soil  pressures 
can  be  resisted  mainly  by  thrust  forces  in  the  conduit  wall.   In  cases 
where  the  soil  support  is  moderate  to  weak  and  in  less  flexible 
conduits,  significant  bending  stresses  accompany  deflection.   Demmin ' s 
tests  showed  tloat  large  bending  stresses  can  develop  after  relatively 
small  deflection  in  long-span  pipe  arches  with  shallow  cover.   Although 
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FIGURE  II   WALL  CRUSHING  APPARATUS  (after  LANE,  1965) 
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FT 

the  pipe  arch  was  very  flexible  (r^  =   1-5),  the  wall  stiffness  (EI) 

R 
was  relatively  large,  and  only  small  changes  in  the  radius  of 

curvature  were  needed  to  produce  large  moments  and  local  yielding. 

The  effects  of  local  yielding  on  conduit  performance  and  its 

interaction  with  the  buckling  mechanism  has  not  been  addressed 

adequately  and  is  presently  a  source  of  major  uncertainty  in  adapting 

experience  with  conduits  of  normal  size  (up  to  12  feet  in  diameter) 

to  the  design  of  long-span  conduits  (up  to  60  feet  or  more  in  span). 

Inelastic  behavior  is  also  important  in  smaller  conduits,  as 

reflected  by  Howard's  (1971)  tests.   In  the  weaker  clay  backfill  the 

most  flexible  conduit  failed  catastrophically  by  elastic  buckling. 

Less  flexible  conduits  exhibited  considerable  inelastic  behavior 

before  collapse.   Initial  yielding  (crimping)  occurred  locally  in 

several  locations  at  approximately  one  half  the  load  capacity. 

Howard  believed  the  crimps  resulted  from  plastic  hinging,  but  because 

deflections  before  crimping  occurred  increased  as  the  conduits 

became  less  flexible,  the  writers  believe  that  local  inelastic 

buckling  initiated  wall  disturbance.   Crimping  occurred  at  larger 

deflections  because  of  increasing  resistance  to  buckling.   Analysis 

of  inelastic  buckling  and  definition  of  subsequent  performance  is 

not  possible  at  this  time,  and  more  study  is  needed  in  this  area.   In 

denser  clay  backfill  initial  crimping  resulted  from  elastic  local 

buckling.   Bending  restraint  provided  by  the  stronger  backfill  induced 

elastic  buckling  at  overpressures  larger  than  those  producing 

inelastic  crimping  in  an  identical  conduit  supported  by  weaker  backfill. 

It  is  not  possible  to  distinguish  visually  between  the  two  buckling 

modes;  therefore,  tests  on  buried  conduits  should  include 
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instrumentation  designed  to  monitor  strain  distributions  before  and 
after  buckling  is  initiated. 

Performance  Limits 

Associated  with  the  difficulties  involved  in  identifying  the 
mode  of  failure  is  the  question  of  whether  or  not  a  given  mode  of 
failure  should  be  considered  as  a  performance  limit,  i.e.  the  load 
at  which  serviceability  of  the  conduit  is  no  longer  acceptable. 
Referring  to  Fig.  12;  it  is  seen  that: 

(1)  a  local  buckle  seldom  indicates  incipient  collapse, 

(2)  a  snap-through  buckle  may,  or  may  not,  indicate  incipient 
collapse,  and 

(3)  wall  crushing  and  seam  distress  can  occur  locally  at 
loads  much  below  the  collapse  load. 

Thus,  it  is  far  from  clear  what  conditions  should  be  considered  as 
performance  limits  for  design  purposes. 

Although  when  local  buckling  or  wall  crushing  has  occurred  the 
conduit  may  not  be  near  collapse,  its  structural  load  capacity  has 
been  impaired.   Over  long  periods  of  time  the  conduit  may  continue 
to  deform  with  further  deterioration  of  its  structural  capacity, 
and  the  possibility  of  damage  from  corrosion  is  enhanced.   In  the 
writers'  view,  studies  are  needed  to  determine  when  local  buckling 
should  be  treated  as  a  performance  limit.   Wall  crushing  and  seam 
separations  involving  breaking  of  welds,  shearing  of  bolt  or  rivet 
holes,  and  the  like,  should  probably  be  considered  as  performance 
limits.   On  the  other  hand,  a  seam  that  is  designed  to  slip  at  a 
predetermined  load,  but  maintain  its  resistance  to  bending,  can 
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serve  to  reduce  and  redistribute  interface  pressures,  and  thereby  enhance 
performance  considerably. 

Role  of  Controlled  Experiments 

The  ability  to  predict  the  critical  loads  causing  elastic 
buckling  or  plastic  yielding  (including  the  development  of  hinges) 
is  still  rudimentary,  and  the  writers  are  not  aware  of  any  methods 
for  predicting  inelastic  buckling  in  buried  conduits.   Thus,  the 
only  way  to  insure  that  a  buried  conduit  will  perform  satisfactorily 
is  to  design  the  conduit  so  that  wall  crushing  will  be  the  critical 
failure  mode.   To  do  this  with  assurance,  the  use  of  well -compacted 
granular  backfill  around  conduits  is  presently  required,  and  requisite 
experience  for  selecting  appropriate  wall  sections  is  available  only  for 
corrugated  steel  culverts  up  to  diameters  of  about  12  feet.   To 
extend  this  experience  to  conduits  of  other  materials  and  larger 
spans,  and  to  other  types  of  backfill,  will  require  accumulation  of 
extensive  performance  records  over  long  periods  of  time.   Alternately, 
prediction  capabilities  can  be  improved  and  verified  using  a  limited 
number  of  controlled  experiments:  the  results  can  then  be  applied  to 
existing  performance  records,  and  improved  design  methods  developed 
in  less  time,  and  at  a  much  lower  cost. 
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V .   SUMMARY 

Despite  the  complexity  of  interaction  between  flexible  conduits 
and  backfill  soil,  and  the  associated  difficulties  in  predicting 
conduit  behavior,  the  performance  record  of  typical  corrugated  steel 
culverts  (less  than  12  feet  diameter)  is  very  good.   The  majority  of 
failures  that  have  occurred  can  be  traced  to  excessive  live  loading 
during  construction  or  to  sub-standard  backfill  soil.   Inability  to 
model  soil  and  conduit  properties  accurately  and  to  identify  the 
conditions  under  which  different  failure  modes  are  critical  has  been 
compensated  for  by  extensive  experience  and  the  use  oC  appropriate 
safety  factors  in  design.   On  the  other  hand,  for  conduit  materials 
other  than  steel,  for  cohesive  backfills,  and  for  longer  spans, 
comparable  empirical  information  is  not  available.   Only  better 
understanding  of  the  possible  failure  mechanisms  can  lead,  in  a  short 
period  of  time,  to  improved  predictive  capabilities  and  hence  to 
more  economical  designs. 

This  report  has  reviewed  many  of  the  theoretical  and  experimental 
studies  that  have  been  conducted  on  buried  flexible  conduits  and  has 
attempted  to  assess  the  value  of  these  studies  for  the  purpose  of 
establishing  rational  design  criteria  for  buried  flexible  conduits. 
The  overall  results  are  summarized  below. 

A.   Wall  Crushing 
Wall  crushing  (yielding  of  the  conduit  section  due  to 
thrust  forces  only)  has  been  observed  experimentally  in  conventional- 
sized  flexible  conduits  supported  by  very  dense  granular  backfill. 
Deflections  before  failure  are  consequently  very  small  which,  coupled 
with  the  relatively  low  stiffness  of  the  section,  insures  that  bending  will 
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contribute  little  to  the  wall  stresses.   If  failure  is  by  wall  crushing, 
then  the  conduit  has  supported  its  overburden  load  in  the  most  efficient 
manner.   The  goal  is  to  select  the  geometry,  culvert  section,  backfill, 
and  construction  methods  to  approach  this  ideal  as  closely  as  possible. 

B.   Seam  Failure 
Seam  separation  is  probably  the  simplest  failure  mechanism  to 
identify.   Seam  integrity  is  a  function  of  the  type  of  culvert,  the 
construction  technique,  and  the  thrust  load.   Strength  tests  by  Lane 
(1965)  showed  that  seam  failure  (either  rivet  shearing  or  weld  rupture) 
is  critical  in  corrugated  culverts  less  than  2  feet  in  diameter.   However, 
tests  on  culverts  with  larger  diameters  (up  to  4  feet)  showed  that  wall 
crushing  is  more  critical.   Watkins  and  Moser's  tests  also  showed  that 
some  other  mode  of  wall  disturbance  preceeds  seam  distress.   Although 
seam  slippage  occurs  at  loads  less  than  the  seam  capacity,  structural 
integrity  of  the  culvert  is  not  immediately  affected  and,  generally, 
a  more  favorable  distribution  of  interface  soil  pressure  results. 
Research  on  seam  design  could  be  very  rewarding:   a  seam  that  can 
continue  to  slip  without  impairing  bending  resistance  at  thrust  loads 
below  those  required  to  induce  wall  crushing  could  result  in  signifi- 
cantly more  economic  designs. 

C.   Elastic  Buckling 
Elastic  buckling,  the  development  of  an  instability  before  yield- 
ing is  initiated  in  the  conduit  section,  can  occur  in  high  modes  if  the 
interface  pressures  are  reasonably  liniform  and  imperfections  or  local 
residual  stresses  are  insignificant.   Otherwise,  elastic  buckling  is 
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initiated  as  a  local  buckle  in  one  (or  a  few)  locations  and  may  occur 
at  the  crown,  invert  or  other  positions  depending  on  the  location  of 
the  critical  combination  of  thrust,     bending  moment,  imperfection  and 
residual  stress.   A  snap-through  buckle  may  be  an  initial  sign  of  dis- 
tress or  it  may  develop  at  higher  loads  after  yielding  or  smaller 
local  crimping  has  occurred  at  another  location.   The  many  high  mode 
buckling  theories  considering  elastic  soil  support  are  relatively 
similar,  excluding  the  different  approaches  for  formulating  the  modulus 
of  soil  support.   Most  expressions  for  the  critical  buckling  pressure 
are  lower-bound  solutions.   However,  even  lower-bound  solutions  generally 
overestimate  the  elastic  buckling  loads  observed  in  controlled  experi- 
ments.  Imperfections  and  residual  stresses  can  lead  to  early  local 
buckling.   Several  methods  have  been  proposed  to  account  for  the  accom- 
panying reduction  in  critical  pressure,  but  none  of  these  methods  has 
been  verified  by  controlled  experiments  of  sufficient  scope  and  generality. 
Of  the  snap-through  buckling  theories,  only  Kloppel  and  Clock  treat  non- 
uniform boundary  pressures  and  consider  the  deflected  conduit  shape. 
Their  solution  is  probably  the  most  satisfactory  representation  currently 
available,  yet  it  may  be  overconservative  in  some  situations  and  unsafe 
in  others.   Materials  with  high  yield  points  offer  advantages  in  over- 
coming the  effects  of  imperfections  and  residual  stresses  on  the  load 
capacity.   Studies  are  needed  to  quantify  the  improvements  that  can  be 
realized  economically. 

Kloppel  and  Clock's  elastic  snap-through  buckling  theory  indicates 
that  elliptical  conduits  are  more  susceptible  to  buckling  than  those 
of  circular  shape.   For  long-span  conduit  and  backfill  stiffnesses  in 
current  use,  the  critical  buckling  load  would  be  approximately  half 
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that  of  an  equivalent  circular  shape.   Thus,  buckling  criteria  adapted 
from  experience  with  circular  culverts  may  not  be  applicable  to  long- 
span  arches  (Stetkar,  1978). 

Critical  buckling  pressures  for  larger  conduits  cannot  be  scaled 
directly  from  small-scale  experiments.   Larger-scale,  controlled  tests 
show  that  elastic  local  buckling  can  be  the  first  sign  of  distress  in 
flexible  conduits  with  moderately  good  soil  support  and  can  occur 
after  relatively  small  deflection.   In  some  cases  the  development  of 
a  local  buckle  constitutes  a  performance  limit,  whereas  in  others  the 
conduit  can  sustain  considerable  additional  load  before  failure  occurs 
(little  resistance  to  further  deflection).   Studies  are  badly  needed 
to  define  the  conditions  under  which  local  buckling  should  be  con- 
sidered as  a  performance  limit . 

D.   Excessive  Deflection 
Excessive  deflections  have  been  observed  in  flexible  conduits  of 
conventional  size  supported  by  weak  soil,  and  in  long-span  conduits 
supported  by  relatively  strong  soil.   In  the  former  case  deflections 
may  exceed  15%  of  the  diameter  before  collapse  is  imminent,  while  in 
the  latter  case  the  deflections  may  be  less  that  5%  of  the  span. 
Failure  by  excessive  deflection  in  long-span  conduits  with  good  backfill 
can  be  triggered  by  bending  stress  concentrations  that  are  produced 
during  construction.   Zones  of  plastic  yielding  can  develop,  and  be- 
cause soil  support  is  low  in  the  vicinity  of  the  crown  (due  to  shallow 
cover) ,  catastrophic  collapse  in  the  form  of  a  snap-through  buckle 
can  follow.   Analytical  models  predicting  the  bending  moment  distribu- 
tion during  construction  of  long-span  conduits  are,  at  best,  rudimentary. 
Thus,  careful  monitoring  of  field  installations  are  needed  to  advance 
the  state-of-the-art. 
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Excessive  deflection  of  typical  culverts  with  deep  burial  does  not 
initially  involve  bending  stress  concentrations  near  the  crown.   De- 
flections at  the  crown  induce  soil  arching,  which  tends  to  equalize 
the  boundary  pressure  distribution.   If  the  required  deflection  is 
sufficiently  large,  flattening  of  the  culvert,  coupled  with  the  develop- 
ment of  plastic  yielding  at  the  springlines,  can  induce  an  inversion 
of  curvature  involving  a  large  fraction  of  the  circumference  in  the 
vicinity  of  the  crown.   If  loading  is  continued,  the  culvert  will 
collapse.   Such  excessive  deflection  will  not  occur  in  granular  back- 
fill densified  to  95%  Std.  AASHO  density:   however,  if  a  cohesive  back- 
fill compacted  to  this  specification  is  used,  more  deflection  is  re- 
quired to  mobilize  equivalent  lateral  soil  resistance  which,  coupled 
with  gradual  loss  of  circumferential  support  due  to  time-dependent 
behavior,  can  lead  to  failure.   More  studies  on  the  behavior  of  con- 
duits supported  by  cohesive  backfills  are  needed. 

E.   Plastic  Hinging 
With  less  flexible  conduits,  plastic  yielding  under  combined 
bending  and  thrust  stresses  can  continue  until  a  plastic  hinge  is  formed 
before  excessive  deflection  at  the  crown  occurs;  that  is,  the  rotational 
resistance  at  a  point  in  the  cross-section  becomes  fully  mobilized. 
If  additional  load  is  applied,  redistribution  of  boundary  pressures 
must  occur  and  if  rotational  resistance  becomes  fully  mobilized  at  a 
sufficient  number  of  locations  the  conduit  will  collapse.   Even  if  a 
collapse  mechanism  does  not  develop,  failure  by  excessive  deflections 
or  inelastic  buckling  may  occur  prematurely,  and  the  conduit  will  not 
be  supporting  load  in  its  most  efficient  manner.   Consequently,  bending 


94. 


moments  can  have  important  influences  on  buried  conduit  behavior.   It 
has  been  shown  experimentally  that  plastic  yieldinj;  may  develop  aTter 
small  deflections,  especially  in  large  diameter  culverts  and  pipe  arclies, 
Finite  element  analyses  can  estimate  the  contribution  of  bending  to  the 
total  wall  stress,  but  present  formulations  are  not  well  adapted  to  the 
analysis  of  post-hinging  behavior.   Analytical  techniques  need  to  be 
extended  to  permit  prediction  of  collapse  loads. 

F.   Inelastic  Buckling 
No  theory  has  been  developed  to  predict  inelastic  buckling 
(instability  after  yielding  has  been  initiated)  of  buried  conduits. 
To  date,  buckling  tests  have  lacked  sufficient  instrumentation  even  to 
develop  empirical  insight  into  the  problem.   Inelastic  buckling  appears 
to  be  a  critical  failure  mode  in  flexible  culverts  whenever  the  de- 
flections necessary  to  mobilize  lateral  soil  support  are  sufficient 
to  induce  yielding,  and  in  less  flexible  culverts  subjected  to  non- 
uniform loading.   Thus,  it  is  suspected  that  inelastic  buckling  may  be 
an  important  failure  mode  in  long-span  culverts.   Probably  no  other 
aspect  of  buried  conduit  behavior  is  in  greater  need  of  additional 
study. 

G.   Recommended  Experimental  Studies 
Because  many  aspects  of  conduit-soil  interaction  cannot  presently 
be  modeled  theoretically,  extensively  instrumented,  controlled  tests 
on  large  scale  conduits  are  badly  needed.   These  tests  should  be 
conducted  with  a  sufficient  range  of  conduit  flexibility  and  soil 
support  (granular  and  cohesive)  to  encompass  all  types  of  failure 
modes.   Loading  should  be  continued  to  collapse  so  that  the 
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relationship  between  initial  wall  disturbance  and  the  collapse  load 
can  be  studied.   In  each  test  the  nature  and  quality  of  soil  support 
should  be  directly  measured  around  the  conduit  and  then  compared  to 
results  obtained  from  triaxial  and  plane-strain  laboratory  tests  at 
comparable  confining  pressures.  , 

Conduits  should  be  instrumented  to  measure  their  deflected  shape, 
the  distribution  of  strains  due  to  thrusts  and  bending  moments,  and 
the  distribution  of  interface  soil  pressures.   This  will  allow  reexam- 
ination of  the  validity  of  simplifying  theoretical  and  design  assumptions, 
and  establish  the  general  applicability  of  the  test  results.   The 
ability  of  finite  element  computer  techniques  to  predict  the  behavior 
of  conduits  under  a  large  variety  of  controlled  conditions  can  also 
be  evaluated.   Measurement  of  shear  stresses  at  the  interface  is  needed, 
as  considerable  doubt  remains  regarding  the  influence  of  traction 
forces  on  conduit  behavior. 

Because  imperfections  and  residual  stresses  can  have  a  large 
effect  on  conduit  behavior,  tests  on  unsupported  conduits  should  be 
conducted  and  compared  with  theoretical  solutions  to  assess  their 
influence  on  the  load  capacity.   Also,  in  tests  where  soil  support  is 
provided,  residual  stresses  due  to  fabrication  should  be  measured  inde- 
pendently so  that  the  true  conduit  stresses  at  failure  can  be  obtained. 
All  these  measurements  will  provide  a  basis  for  further  developments 
in  the  design  of  all  kinds  of  buried  conduits,  including  that  of 
long-span  arches. 
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